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. Abstract: We propose to use the Mt2 concept to measure the masses of all particles in 

1/-^ . SUSY-hke events with two unobservable, identical particles. To this end we generalize the 

' usual notion of Mt2 and define a new Mj^2^'^^ variable, which can be applied to various 

subsystem topologies, as well as the full event topology. We derive analytic formulas for 
(yr^ \ its endpoint M^P^^^^x ^ ^ function of the unknown test mass Mc of the final particle in 

' the subchain and the transverse momentum pr due to radiation from the initial state. We 

show that the endpoint functions Mj?2^^axi^c, Pt) may exhibit three different types of kinks 
^ ' and discuss the origin of each type. We prove that the subsystem M^^''^'^ variables by 

' themselves already yield a sufficient number of measurements for a complete determination 

of the mass spectrum (including the overall mass scale). As an illustration, we consider 
the simple case of a decay chain with up to three heavy particles, X2 Xi — > Xq, which 
is rather problematic for all other mass measurement methods. We propose three different 
MT2-based methods, each of which allows a complete determination of the masses of particles 
Xq, Xi and X2- The first method only uses M^g'^''^^ endpoint measurements at a single fixed 
value of the test mass Mc- In the second method the unknown mass spectrum is fitted to 
one or more endpoint functions M!f!^^^^^{Mc,PT) exhibiting a kink. The third method is 
hybrid, combining Mt2 endpoints with measurements of kinematic edges in invariant mass 
distributions. As a practical application of our methods, we show that the dilepton W~^W~ 
and ti samples at the Tevatron can be used for an independent determination of the masses 
of the top quark, the W boson and the neutrino, without any prior assumptions. 
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1. Introduction 

The ongoing Run II of the Fermilab Tevatron and the now commencing run of the Large 
Hadron Collider (LHC) at CERN are on the hunt for new physics beyond the Standard 
Model (BSM) at the TeV scale. Arguably the most compelling phenomenological evidence 
for BSM particles and interactions at the TeV scale is provided by the dark matter problem [1], 
whose solution requires new particles and interactions BSM. A typical particle dark matter 
candidate does not interact in the detector and can only manifest itself as missing energy. 
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Figure 1: The generic event topology under consideration in this paper. The particles Xi, 1 < i < 
are new BSM particles which appear as promptly decaying, on-shell intermediate resonances. The 
particles Xi are the corresponding SM decay products, which are all visible in the detector, i.e. we 
assume that there are no neutrinos among them. ISR stands for generic initial state radiation with 
total transverse momentum pV- is a BSM particle which is invisible in the detector. The integer 
n counts the total number of intermediate BSM particles in each chain, so that the total number of 
BSM particles in each chain is n + 1. For simplicity, in this paper we shall only consider symmetric 
events, in which the two decay chains are identical. The generalization of our methods to asymmetric 
decay chains is straightforward. 

The dark matter problem therefore greatly motivates the study of missing energy signatures 
at the Tevatron and the LHC [2]. 

The long lifetime of the dark matter particle is typically ensured by some new exact 
symmetry^, under which the SM particles are neutral, while the BSM particles are charged. 
This setup implies that the new particles will be pair-produced, and each of the two cascades 
will terminate in the dark matter candidate, giving rise to missing energy in the detector. (A 
generic example of this topology is shown in Fig. [l|.) Since the energies and momenta of the 
final two invisible particles Xq are not measured, one cannot directly apply resonance mass 
reconstruction techniques^ . This represents a significant challenge for determining the masses 
Mi of the new particles Xi . In recognition of this problem, there has been a recent resurgence 
of interest in the development of different methods for mass measurements in cascade decays 
with missing energy [12-48]. Most of these techniques fall into one of the following three 
categories: 

• I. Endpoint methods. They rely on the kinematic endpoints [12,14-16,20,21,25] or 
shapes [22,23,46] of various invariant mass distributions constructed out of the visible 
(SM) decay products Xi in the cascade chain. 

• II. Polynomial methods. Here one attempts exact event reconstruction using the 
measured momenta of the SM particles and the measured missing transverse momentum 

^Some well known examples are: i?-parity in supersymmetry [3], KK parity in models with extra dimensions 
[4-7], T-parity in Little Higgs models [8,9], f/-parity [10,11] etc. 
■^See, however, Section 2.2 
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[18,19,27,38,39]. 

• III. methods. These methods explore the transverse invariant mass variable Mt2 
originally proposed in [13] and later used and developed in [17,24,26,28,37,40,41,45]. 
Recently it was shown that under certain circumstances, the endpoint of the Mt2 
distribution, when considered as a function of the unknown test mass Mq of the lightest 
new particle Xq, exhibits a kink and the true mass Mq of Xq, i.e. at Mq = Mq [29-32,36]. 

One could also combine two or more of these techniques into a hybrid method, e.g. a mixed 
polynomial and endpoint method [34], a mixed Mt2 and endpoint method [33,43], or a mixed 
Mt2 and polynomial method [47,48]. In Section ^ we shall describe in detail each of these 
three basic approaches I - III. We shall then contrast them to each other and discuss their 
pros and cons. In particular, we shall concentrate on their applicability as a function of the 
length of the decay chain, i.e. the number n of intermediate resonances in Fig. I. We shall 
find that for sufficiently long decay chains, namely n > 3, each method I - III hy itself is able 
to completely determine the unknown particle spectrum, at least as a matter of principle. 
Therefore, if Nature is so kind to us as to present us with such a long decay chain, it does not 
really matter which of the three methods above we decide to use - sooner or later, success 
will be guaranteed with each one. 

However, if the decay chain happens to be relatively short, i.e. n < 2, neither method 
I, nor method II, nor a hybrid combination of I and II will be able to completely determine 
the unknown particle mass spectrum. In contrast, method III by itself can still provide a 
sufficient number of measurements for a complete determination of the mass spectrum of the 
new particles. We argue that in order to achieve this, the conventional Mt2 variable needs to 
be promoted to a more general quantity M^g'^'^^ , which can be applied not only to the whole 
event, but also to a particular sub-chain starting at Xp and ending in Xc [41,45]. We present 
the basic steps for this generalization in Section ^, where we also introduce our conventions 
and notation. Then in Section Q we concentrate on the problematic case of n < 2 and discuss 
what type of M^^'^^ measurements are available in that case. We then show that the newly 
defined M^f^^'^^ may also exhibit a kink in the graph of its endpoint m!^^^^ as a function 
of the test mass M^- In order to be able to properly interpret this kink, we derive analytic 
expressions for the function M^Ym^xi^c ,Pt), including the effect of initial state radiation 
(ISR) with some arbitrary transverse momentum px (see Fig. |l|). In all cases, a kink would 
always appear at Mc = M^. 

g <'ri(^c,PT) \ , ( dM!p-dAMc,PT, 

dM, ) ^1 dM, ' ' ^ ^ 

/ Mc=Mc-e \ / Mc=Mc+€ 

and the value of that point reveals the true mass Mp of the mother particle Xp-. 

K7d.iMc,PT) = Mp . (1.2) 



However, there may be up to three different reasons for the origin of the kink (O). For 
example, in the case of -^^2 max(-^O)PT) and M!jp2'^^^^{Mi,px) with non-zero px, the kink 
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arises due to recoils against the ISR jets, as explained in [30, 31] (see Sections 4A and O 
below). On the other hand, in the case of Mj,^ max(-^o,PT = 0), the kink is due to the variable 
mass of the composite system of SM particles {X1X2}, as already observed in Refs. [29,32] 
(see Section below). Finally, for Mg^;"]^ (Mq , = 0), we encounter a new type of kink, 
which arises due to the decay of a heavier particle (in this case X2) upstream (see Section 



4J). 

In Section [s] we propose three different methods for measuring the masses of all the 
particles in the problematic case of an n = 2 decay chain. With the first method, presented 
in Section |5.1| , we always consider a fixed value of the test mass (for convenience we choose it 



to be zero), and perform a sufficient number of M^f^^'^^ endpoint measurements for various n. 



p and c. In our second method, described in Section |5.2| , we choose a suitable M^*^'^'^^ variable 



whose endpoint M^P^^^ax exhibits a kink, and then fit for the function -^j^'max(-^c)PT)- Our 
last method, presented in Section [5^ , is hybrid, in the sense that we use combined information 
from the measured endpoints of some Mt2 distributions, as well as the measured endpoints 
of certain invariant mass distributions. Neither of our three methods relies on reconstructing 
the actual momentum of each missing particle. 

All of our discussion throughout the paper will be completely model-independent and 
can be applied to any BSM scenario, including super symmetry, extra dimensions, little Higgs 
theory etc. In Section |5|, however, we shall use a specific example in order to illustrate each 
of our three proposed methods. Instead of considering a decay chain of some BSM model, 
we chose to select an example which is already present in the Tevatron data, and will soon 
be tested at the LHC as well: the dilepton event samples from top quark pair production 
and from VF-pair production. Those two dilepton samples satisfy all of our assumptions, and 
would be a perfect testing ground for any new ideas about mass measurements in missing 
energy events from new physics. In Section ^ we will show that using any one of our three 
-MT2-based mass measurement methods, one can in principle determine the mass of each of the 
three particles: top quark, T^-boson, and neutrino, independently and in a completely model- 
independent fashion. Section R contains a summary and a discussion of our main results. In 
Appendix]^ we collect all relevant formulas for the endpoint functions MT2^max^^<^^PT)- 

2. Mass measurement methods in missing energy events 

Let us now discuss in some detail each of the three basic methods I - III for mass measure- 
ments in missing energy events. The basic topology is shown in Fig. |l|, where the particles 
< i < n, (denoted in red) are new BSM particles, and the particles Xi, 1 < i < n, 
(denoted in black) are the corresponding SM decay products. ISR stands for generic initial 
state radiation with total transverse momentum pT- Xq is a dark matter candidate which 
is invisible in the detector. For simplicity, in this paper we shall make two assumptions, 
each of which can be easily relaxed without significantly changing our conclusions. First, we 
shall assume that the intermediate particles Xi, 1 < i < n, are all on-shell, i.e. their masses 
Mi = obey the hierarchy 

M„ > Mn-i > . . . > Ml > Mo . (2.1) 



-4- 



Consequently, all decays along our decay chain are two-body, i.e. each SM decay product Xi 
is a single particle of mass rrii = nixi- ^^^^ paper we will also only concentrate on the 
commonly encountered case where Xi is either a lepton, photon or jet, i.e. massless: 

mi = 0, i = l,...,n. (2.2) 

Second, we shall also assume that our events are symmetric, i.e. the two decay chains are 
identical. Again, this assumption can be easily dropped and one could consider asymmetric 
events as well (see, for example, Ref. [45]). These assumptions are made only for simplicity. 
All of our subsequent discussion can be easily generalized to include off-shell decays, simply 
by promoting some of the visible SM particles Xi to composite particles with variable mass. 
For example. Appendix ^ already contains some results for an off-shell case, while the study 
of asymmetric events is postponed for future work [49]. 

We shall use the integer n to count the total number of intermediate on-shell BSM 
particles in each chain. Then, the total number of new BSM particles in the decay chain is 
n + 1. With those preliminaries, we are ready to discuss each of the three different methods 
for mass measurements in missing energy events. 

2.1 Endpoint method 

With this method, one forms the invariant mass distributions Mx^^xi^.^xi^^ of various groups 
of k SM decay products Xi, where the number k in principle can range from 2 to n. Each 
such distribution exhibits an upper kinematic endpoint, which can be related to the under- 
lying unknown masses Mi. If one makes a sufficient number of independent upper endpoint 
measurements, the system of equations giving the kinematic endpoints Ej in terms of the 
masses Mj 

Ej=Ej{Mo,...,Mn), j = 0,...,n (2.3) 

can be solved for the masses Mj, although on some occasions the solution may not be unique 
- see, e.g. [50, 51]. 

Clearly, the method will be fully successful only if the number of measurements Nm is no 
less than the number of unknown parameters Np. For the decay chain of Fig. |l|, the number 
of unknown mass parameters Np is simply the total number of BSM particles: 

Np = n + 1 . (2.4) 

How many measurements Nfn are available with this method? The answer to this question 
depends on the length of the decay chain. It is easy to see that, if n = 1, there are no 
endpoint measurements at all; if n = 2, there is a single measurement of the endpoint of the 
Mx^x2 distribution, etc. In general, for an arbitrary fixed n, the number of different invariant 
mass distributions Maj-^a;.^..,^,.^ that one can form and study, is equal to the number of ways 
in which we can select a group of at least two objects from a set of n objects, and is given by 

iV^ = 2" - (n + 1) . (2.5) 
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Strictly speaking, eq. ( |2.5D only gives an upper bound on the number of independent upper 
endpoints in the invariant mass distributions. Indeed, there are certain cases where not all 
of the upper kinematic endpoints are independent. For example, consider the familiar case 
of a squark decay chain in supersymmetry: X3 = q, X2 = — ^ ™d Xq = Xi- The SM 

decay products consist of a quark jet q and two leptons: i'^ and £~ . It is well known that in 
some regions of parameter space the upper endpoint of the Mgu distribution does not provide 
an independent measurement, since it can be related to the upper endpoints of the M^+n- 
and Mg£(/jjg/j) = max{Mg£+,Mg£-} distributions [20]. Fortunately, one can use additional 
measurements from the lower endpoints of suitably restricted invariant mass distributions, 
e.g. MqiK^Q^T^-^ [14]. We see that the precise count of the number of measurements Nm available 
in the endpoint method is somewhat model-dependent, but nevertheless, the estimate 
is sufficient to make our main point below. 

From eqs. (^^) and (|2.5| ) it readily follows that the number of undetermined parameters 
with this method is 

Np-N^ = 2(n + 1) - 2" . (2.6) 

The dependence of this quantity on the length of the decay chain is plotted in Fig. § with 
red open circles, connected with red line segments. The yellow-shaded region in the figure 
is where Nm > Np, so that we have a sufficient number of measurements for a complete 
determination of the heavy particle spectrum. Conversely, whenever a symbol appears inside 
the white region, where Nm < Np, there is only partial information about the mass spectrum 
and the spectrum cannot be fully determined. 

Fig. ^ reveals that the endpoint method cannot succeed unless n > 3. This conclusion has 
already been confirmed by numerous studies of various low-energy SUSY models, where one 
considers a decay chain of sufficient length: n = 3 as in the squark example mentioned above, 
or n = 4 as for a gluino chain [21]. On the other hand, if n = 1 or n = 2, with this method we 
are unable to pin down all of the new particle masses, even as a matter of principle. These 
are exactly the cases where the additional information from mass measurements at future 
lepton colliders has been seen as extremely useful [52]. 

2.2 Polynomial method 

The basic idea behind the method is to use all of the available experimental information in 
each event, and enforce a sufficient number of constraints, which would allow to actually solve 
for the unknown momenta of the missing particles Xq. Before we analyze this method in more 
detail, let us introduce some of our notations and conventions. We shall use lowercase letters 
to denote various quantities relating to the SM particles Xi, 1 < i < n. At the same time, 
we shall use a superscript (k) to denote whether a particular quantity belongs to the upper 
{k = 1) or lower {k = 2) decay chain in Fig. Il|. For example, the 4-momentum of the SM 

(k) 

particle Xi in the k-th. chain will be denoted as p\ , the corresponding transverse momentum 

—(k) 

will be Pj^rp , while the mass of Xi will simply be m^. On the other hand, uppercase letters 
will denote quantities relating to the BSM particles Xi, < i < n. Thus the 4-momentum 

(k) 

of the BSM particle Xi appearing in the k-th chain is , the corresponding transverse 
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Figure 2: The dependence of the number of undetermined parameters Np — Nm as a function 
of the number n of intermediate heavy resonances in the decay chains of Fig. ^ for various mass 
determination methods: Mt2 method (green, open squares), endpoint method (red, open circles), 
polynomial method for N^v — 2 (blue, x symbols), or a hybrid method which is a combination of the 
latter two methods (magenta, symbols). Within the yellow-shaded region the number of unknowns 
Np does not exceed the number of measurements Nm for the corresponding method, and the mass 
spectrum can be completely determined. 

~'(k) 

momentum is P^rp , and the mass is Mj. One should keep in mind that for SM particles the 
index i runs from 1 to n, while for BSM particles i runs from to n. 

In order to apply the polynomial method, one uses the experimentally measured 4- 

(k) -, 

momenta pi as well as the missing transverse momentum pT,miss in the event. Then, one 
imposes the mass shell constraints for the intermediate BSM particles Xi and tries to solve 

(k) 

the resulting system of equations for the 8 unknown components of the 4-momenta Pq of 
the missing particles Xq. Including the n + l unknown masses Mj, this amounts to a total of 



unknown parameters. How many measurements (constraints) are present in this case? First, 
there is a total of 2(n + 1) mass-shell conditions: one for each BSM particle Xi in each of the 




(2.7) 



two decay chains in Fig. |T] 




i = 0, 1, . . . , n . 




Using energy and momentum conservation 



pW^pW + ^^W^ A. = 1,2, 



(2.9) 



(k) 

these constraints can be rewritten in terms of the unknown variables Pq : 

m! = (^pI,'^ +±pf^ =(^^o'^+&rj ' ^ = 0,l,...,n. (2.10) 

Furthermore, the measurement of the missing transverse momentum pT,miss provides two 
additional constraints 

-^OT + -^OT ~ PT,miss (2-11) 
~'{k) 

on the unknown transverse momentum components P^q^ ■ Therefore, the total number of 
measurements is 

Ar^ = 2(n + l) + 2 = 2n + 4 (2.12) 

and the number of undetermined parameters for any given event is readily obtained from 
^ and ( pl2D 

Np- Nm = fy-n . (2.13) 

However, one might do better than this, by combining the information from two or more 
events [18, 19]. For example, consider another event of the same type. Since the n + 1 
unknown masses were already counted in eq. (^]^), the second event introduces only 8 new 
parameters (the 4- momenta of the two Xq particles in the second event), bringing up the 
total number of unknowns in the two events to 

ATp = 8 + 8 + (n + 1) = n + 17 . (2.14) 

At the same time, the constraints ( p.lOj ) and ( |2.11| ) are still valid for the second event, which 
results in 2n + 4 additional constraints. This brings the total number of constraints to 

iV„ = (2n + 4) + (2n + 4) = 4n + 8 . (2.15) 

Subtracting ( |2.14| ) and ( p. 151 ), we get 

iVp - iV„ = 9 - 3n . (2.16) 

Comparing our previous result (|2ll ) with ( ^lel) , we see that the latter decreases much 
faster with n, therefore, when using the polynomial method, combining information from two 
different events is beneficial for large n (in this example, for n > 3). 

Following the same logic, one can generalize this parameter counting to the case where 
the polynomial method is applied for sets of N^v different events at a time. The number of 
unknown parameters is 

Np = n + l + 8Nev , (2.17) 

the number of constraints is 

= (2n + A)Nev , (2.18) 
and therefore, the number of undetermined parameters is given by 

Np- Nm = n+l-2{n- 2)Nev ■ (2.19) 
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For Nev = 1 and N^v = 2 this equation reduces to ( p. 13 ) and ( 2.16| ), respectively. What is 



the optimal number of events N^v for the polynomial method? The answer can be readily 
obtained from eq. ( 2.ig| ), where N^v enters the last term on the right-hand side. If this term 



is negative, increasing N^.^ would decrease the number of undetermined parameters, therefore 
it would be beneficial to combine information from more and more different events. From 
eq. (2.19) we see that this would be case if the decay chain is sufficiently long, i.e. n > 3. 



On the other hand, when n = 1, considering more than one event at a time is actually 
detrimental - we are adding more unknowns than constraints. In the case of n = 2, the 
number of undetermined parameters Np — Nm is actually independent of Nev and one might 
as well consider the simplest case of Nev = 1- 

Let us now analyze how successful the polynomial method will be for different decay chain 
lengths. The number of undetermined parameters (|2.16| ) for N^v = 2 is plotted in Fig. ^ with 
blue X symbols, connected with blue line segments. We see that the polynomial method will 
be successful in determining all the masses of the BSM particles only if n > 3. For n = 1 or 
n = 2, there will not be enough measurements for a complete mass determination'^ and the 
best one can do in that case is to obtain a range o/ possible values for the masses Mq, Mi and 
M2 [27]. Recall that in the previous subsection we reached a similar conclusion regarding the 
endpoint method. Therefore, we see that both the endpoint and the polynomial methods, 
when used in isolation, would fail whenever the decay chain is rather short: n = 1 or n = 2. 
This represents a definite problem, since there is no guarantee that the new physics would 
exhibit a long (n > 3) decay chain. Therefore it is worth investigating whether there is an 
alternative method which would be successful in those two cases, i.e. n <2. 

One immediate idea which comes to mind is to use a hybrid method, i.e. combining the 
techniques of the polynomial and endpoint methods [34]. The parameter count in that case 
is very easy to do. The number of unknown parameters is the same as in the polynomial 
method: 

Np = n + l + 8Nev . (2.20) 



Now, however, we need to account for the extra measurements (2.5) which are available from 
the endpoint method. Therefore, the total number of measurements for a hybrid method of 
this type is the sum of (^) and (|2.18|): 



iV^ = 2" - (n + 1) + (2n + 4)Nev ■ (2.21) 
Subtracting (pl20|) and ( ^l2l|) , we get 



Np-Nm = 2{n + 1) - 2" - 2{n - 2)Nev ■ (2.22) 

In Fig. this quantity is plotted for fixed N^v = 2 with magenta (X" symbols. We see that, 
even though the hybrid method performs better than the individual endpoint and polynomial 
methods, it still cannot solve the problem of masses for n = 2\ Therefore, a different approach 



^As can be seen from the more general expression ( ^.19| ), this conclusion will not change even if we consider 
arbitrarily large number of events Ne.v 
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is needed. We shall now argue that the Mt2 method might just provide the solution to this 
problem in the n = 2 case. What is more, we shall show that the Mt2 method can do that all 
hy itself, without using any information derived from the endpoint or polynomial methods. 

2.3 Mt2 method 

Here the number of unknown parameters is still 

Np = n + l . (2.23) 

In the next section we shall prove that, once we consider the notion of subsystem M^2^''^\ 
the total number of MT2-type endpoint measurements is 

Nm = \n{n + l){n + 2) . (2.24) 

D 

Then the number of undetermined parameters with this method is 

Np-N,n = \{n+l)(Q-2n-r?) , (2.25) 

which is plotted with green square symbols in Fig. |2[ We see that for n < 3, the Mt2 method 
is by far the most powerful, and more importantly, it is the only method which is able to 
handle the problematic case of n = 2! 

3. Defining a subsystem Mt2 variable 

The idea for a subsystem Mt2 was first discussed in [41] and applied in [45] for a specific 
supersymmetry example (associated squark-gluino production and decay). Here we shall 
generalize that concept for a completely general decay chain. For this purpose, let us redraw 
Fig. H as shown in Fig. ^. The subsystem Mt2 variable will be defined for the subchain inside 
the blue (yellow-shaded) box in Fig. |. Before we give a formal definition of the subsystem 
Mt2 variables, let us first introduce some terminology for the BSM particles appearing in the 
decay chain. We shall find it convenient to distinguish the following types of BSM particles: 

• "Grandparents" . Those are the two BSM particles Xn at the very top of the decay chains 
in Fig. ^. Since we have assumed symmetric events, the two grandparents in each event 
are identical, and carry the same index n. Of course, one may relax this assumption, 
and consider asymmetric events, as was done in [36,45]. Then, the two "grandparents" 
will be different, and one would simply need to keep track of two separate grandparent 
indices n^^) and n^^). 

• "Parents". Those are the two BSM particles Xp at the top of the subchain used to 
define the subsystem Mt2 variable. In Fig. ^ this subchain is identified by the blue 
(yellow-shaded) rectangular box. The idea behind the subsystem Mt2 is simply to 
apply the usual Mt2 definition for the subchain inside this box. Notice that the Mt2 
concept usually requires the parents to be identical, therefore here we will characterize 
them by a single "parent" index p. 
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Figure 3: An alternative representation of Fig. which iUustratcs the meaning of the subsystem 
variable defined in eq. (|3.3|). 



• "Children". Those are the two BSM particles Xc at the very end of the subchain used to 
define the subsystem Mt2 variable, as indicated by the blue (yellow-shaded) rectangular 
box in Fig. ^. The children are also characterized by a single index c. In general, the 
true mass Mc of the two children is unknown. As usual, when calculating the value of 
the Mt2 variable, one needs to choose a child "test" mass, which we shall denote with 
a tilde, Mc, in order to distinguish it from the true mass Mc of Xc- 

• Dark matter candidates. Those are the two stable neutral particles Xq appearing at 
the very end of the cascade chain. We see that while those are the particles responsible 
for the measured missing momentum in the event (see eq. ( 2.1l| )), they are relevant for 
Mt2 only in the special case of c = 0. 

With those definitions, we are now ready to generalize the conventional M'r2 definition 
[13, 17]. From Fig. | we see that any subchain is specified by the parent index p and the 
child index c, while the total length of the whole chain (and thus the type of event) is given 
by the grandparent index n. Therefore, the subsystem Mt2 variable will have to carry those 
three indices as well, and we shall use the notation M^^''^\ In the following we shall refer to 
this generalized quantity as either "subsystem" or "subchain" Mt2- It is clear that the set 
of three indices {n,p,c) must be ordered as follows: 



n > p > c> . 



(3.1) 



We shall now give a formal definition of the quantity My2'^'^^ j generalizing the original idea of 
Mt2 [13, 17]. The parent and child indices p and c uniquely define a subchain, within which 
one can form the transverse masses M^^^ and M^?^ of the two parents: 



IVlrp yPp , Pp_-^^, . . . , P^^^, I-^rp , IVlc), 



Ak) p(fc). 



= 1,2 . 



(3.2) 



(k) 

Here p\ , c+ 1 < i < p, are the measured 4-momenta of the SM particles within the subchain, 
pj^!^ are the unknown transverse momenta of the children, while Mc is their unknown (test) 
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mass. Then, the subsystem M!j!^^''^^ is defined by minimizing the larger of the two transverse 
masses ( |3.2| ) over the allowed values of the children's transverse momenta Pj:^ : 

M^f'''\Mc)= min |max|Mj!\M|^H| , (3.3) 



where pT indicates any additional transverse momentum due to initial state radiation (ISR) 
(see Figs, pi and ISl). Notice that in this definition, the dependence on the grandparent index 
n enters only through the restriction on the children's transverse momenta P^rp ■ Using 
momentum conservation in the transverse plane 



2 2 n 

fc=l k=l j=l 

~'(k) 

we can rewrite the restriction on the children's transverse momenta as 

E = E 4? + E EpIt = PT,rmss + E E^^^ ' (3-5) 

k=l k=l k=lj=l k=lj=l 



where in the last step we used eq. ( p.ll ). Eq. ( |3.5| ) allows us to rewrite the subsystem Mj!^'^'"^^ 



definition (|3.3D in a form which does not manifestly depend on the grandparent index n: 



{n,p,c) 



T2 



(^^^) = ^. . . ^^^^{max{M«,Mf}} . (3.6) 



V-^ p(f^)-^^ . I Y-2 v-c -.(fc) 
Z^fc=l ^cT —PT,rmss+2^k = l 2^j = l PjT 



However, the grandparent index n is still implicitly present through the global quantity 
PT,miss, which knows about the whole event. We shall see below that the interpretation of 
the experimentally observable endpoints, kinks, etc., for the so defined subsystem Mpg'^'"^^ 
quantity, does depend on the grandparent index n, which justifies our notation. 

We are now in a position to compare our subsystem M!j!^^'^^ quantity to the conventional 
Mt2 variable. The latter is nothing but the special case of n = p and c = 0: 

Mt2 = m|"2"'°) , (3.7) 

i.e. the conventional Mt2 is simply characterized by a single integer n, which indicates the 
length of the decay chain. We see that we are generalizing the conventional Mt2 variable in 
two different aspects: first, we are allowing the parents Xp to be different from the particles 
Xn originally produced in the event (the grandparents), and second, we are allowing the 
children Xc to be different from the dark matter particles Xq appearing at the end of the 
cascade chain and responsible for the missing energy. The benefits of this generalization will 
become apparent in the next section, where we shall discuss the available measurements from 
the different subsystem m!^!^^'^^ variables. 



In conclusion of this section, let us derive the result ( |2.24 ) used in Section 2.3. We count 



how many different subsystem M^'^'^''^^ quantities ( |3.3| ) exist for a given maximum value n 
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of the grandparent index. First, pick a parent index p, which can range from 1 to n. Then, 
for this fixed value of p, the child index c can take a total of p values: < c < p — 1, while 
the grandparent index can take^ a total of n — p + 1 values. Therefore, the total number of 
allowed combinations (n,p,c) is 

n ^ 

yp(n-p+l) = -n(n + l)(n + 2) , (3.8) 

p=l 

in agreement with ( p. 24 ). 



4. A short decay chain X2 ^ Xi ^ Xq 

As we already discussed in Sec. ^, a relatively long (n > 3) new physics decay chain can be 
handled by a variety of mass measurement methods, and in principle a complete determi- 
nation of the mass spectrum in that case is possible at a hadron collider. We also showed 
that a relatively short (n = 1 or n = 2) decay chain would present a major challenge, and a 
complete mass determination might be possible only through Mt2 methods. From now on 
we shall therefore concentrate only on this most problematic case of n < 2. 

First let us summarize what types of subsystem M!f^^'^^ measurements are available in 
the case of n < 2. According to eq. ( |3.8D , there exist a total of 4 different Adff^^'^'^ quantities, 
which are illustrated in Fig. ^. Each M!j?2^''^^ distribution would exhibit an upper endpoint 
-^T^'maa;' "^^osc measurement would provide one constraint on the physical masses. In order 
to be able to invert and solve for the masses of the new particles in terms of the measured 
endpoints, we need to know the analytical expressions relating the endpoints ^^2'^"^^^ to the 

physical masses Mj. In this section we summarize those relations for each M^g'^'^^ quantity 



with n < 2. Some of these results (e.g. portions of Sees. 4.1 and Sees. O) have already 



appeared in the literature, and we include them here for completeness. The discussion in 



Sees. and Sees. on the other hand, is new. In all cases, we shall allow for the presence 
of an arbitrary transverse momentum px due to ISR. This represents a generalization of all 
existing results in the literature, which have been derived in the two special cases px = [32] 
or pt = 00 [31]. 

We shall find it convenient to write the formulas for the endpoints M max ™ terms 
of the actual masses, but in terms of the mass parameters 

Mn f , , 



''Note that different values of the grandparent index n correspond to different types of events. For example, 
in order to study the Myj^'^''^' variables for some given ni, we must look at events of pair-production, 
while in order to form the M^^''^''^^ distributions for another value n2 < ni, we must look at events of X„2 



pair-production. Because of the mass hierarchy (2J), the observation of events of the former type (^ni pair- 
production) guarantees that the collider will eventually be able to also produce events of the latter type {Xn^ 
pair-production). Therefore, the relevant integer for our count is the maximum value of n achievable at a 
given collider experiment. 
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The advantage of using this shorthand notation wih become apparent very shortly. Notice 
that not all of the fi parameters defined in (W^ are independent. For a given maximum 



value of n, the total number of /i parameters from (4.1) is the same as the total number of 



subsystem Mt2 variables and is given by (2.24). All of those fi parameters are functions of 



just n + 1 masses Mj, < i < n, as indicated by eq. ( |2.23| ). Therefore, the fi parameters must 
obey certain relations, whose number is given by ( 2.25| ). For example, for n < 2, we have a 



total of four fi parameters: ^(i^ifi), Ai(2,i,o)i /^(2,2,o) ^iid /i{2,2,i)) and only three masses: Mq, 
Ml and M2, so that there is one constraint: 

/^(2,i,o) (^(2,2,0) - ^(2,2,1)) = t-i-fixo) ■ ('^•2) 
4.1 The subsystem variable M^g^ '^^ 

We start with the simplest case of n = 1 shown in Fig. ^(a). Here ^^^2^'^^^ is the only possibil- 
ity, and it coincides with the conventional Mt2 variable, as indicated by (^]^). Therefore, the 
previous results in the literature which have been derived for the conventional Mt2 variable 
(|3.7| ), would still apply. In particular, in the limit of pr = 0, the upper endpoint M^gl^ix 
depends on the test mass Mq as follows [32] 



^T2il(^o,PT = 0) = ^(1,1,0) + V/^(i,i,o)+^o , (4.3) 



where the parameter /i(i^i^o) is defined in terms of the physical masses Mi and Mq according 
to eq. (|t|) 

Ml) ~ 2Mi 



_Mi ( Mi\_ Ml-Mi 
^(1,1,0) = ^ 1 - 172 - nn. ■ (4-4) 



As usual, the endpoint (4.3) can be interpreted as the mass Mi of the parent particle Xi, so 



that eq. (4.3) provides a relation between the masses of Xq and Xi. In the early literature 



on Mt2, this relation had to be derived numerically, by building the Mt2 distributions for 
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different values of the test mass Mq, and reading off their endpoints. Nowadays, with the 
work of Ref. [32], the relation is known analytically, and, as seen from (4.3), is parameterized 



by a single parameter //(i^i^o)- Therefore, in order to extract the value of this parameter, we 
only need to perform a single measurement, i.e. we only need to study the Mt2 distribution 
for one particular choice of the test mass Mq. We shall find it convenient to choose Mq = 0, 
in which case eqs. (4.3) and ( |4.4] ) give 



<2il(^o = 0,pT = 0) = 2^(,,i,o) = , (4.5) 

providing the required measurement of the parameter /i(i^i^o)- Eq. ( [4.5]) demonstrates the 
usefulness of the M'T2 concept - just a single measurement of the endpoint of the Mt2 
distribution for a single fixed value of the test mass Mq is sufficient to provide us with one 
constraint among the unknown masses (Mi and Mq in this case). 

Unfortunately, one single measurement (^]^) is not enough to pin down two different 
masses. In order to measure both Mq and Mi, without any theoretical assumptions or preju- 
dice, we obviously need additional experimental input. From the general expression (|4.3|) it 



is clear that measuring other m!ix endpoints, for different values of the test mass Mq, will 
not help, since we will simply be measuring the same combination of masses i g) over and 
over again, obtaining no new information. Another possibility might be to consider events 
with the next longest decay chain (n = 2), which, as advertised in the Introduction and 
shown below in Section will be able to provide enough information for a complete mass 
determination of all particles Xq, Xi and X2- However, the existence and the observation 
of the n = 2 decay chain is certainly not guaranteed - to begin with, the particles X2 may 
not exist, or they may have too low cross-sections. It is therefore of particular importance to 
ask the question whether the n = 1 process in Fig. ^(a) alone can allow a determination of 
both Mq and Mi. As shown in Ref. [31], the answer to this question, at least in principle, is 
"Yes", and what is more, one can achieve this using the very same Mt2 variable M^2^''^\ 

The key is to realize that in reality at any collider, and especially at hadron colliders 
like the Tevatron and the LHC, there will be sizable contributions from initial state radiation 
(ISR) with nonzero pT, where one or more jets are radiated off the initial state, before 
the hard scattering interaction. (In Figs. [l|, ^ and ^ the green ellipse represents the hard 
scattering, while "ISR" stands for a generic ISR jet.). This effect leads to a drastic change 
in the behavior of the -^^2 maa;(-^0'-P'j^) function, which starts to exhibit a kink at the true 
location of the child mass Mq = Mq: 

' dM^'£±iMo^\ I dM^'^l{MQ,pT) 

OMq 1 I dMo ' 

/ Mo=Mo-e \ / Mo=Mo+e 

and furthermore, the value of M,^2^^ix ^^^^ point reveals the true mass of the parent as 
well: 

<2il(^o = Mq,pt) = Ml . (4.7) 
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This kink feature (4.6,4.7) was observed and illustrated in Ref. [31] (see their Sec. 4.4). We 



find that it can also be understood analytically, by generalizing the result (4.3) to account for 
the additional ISR transverse momentum pT- Recall that eq. (4.3) was derived in Ref. [32] 
under the assumption that the missing transverse momentum due to the two escaping particles 
Xq is exactly balanced by the transverse momenta of the two visible particles xi used to form 



M 



(1,1,0). 



T2 



-^OT "I" -^dr ~^ Pit ~^ Pit ~ ^ • 



(2) 



(4.8) 



We may sometimes refer to this situation as a "balanced" momentum configuration^. In the 
presence of ISR with some non-zero transverse momentum pT, eq. ( [4.8| ) in general ceases to 
be valid, and is modified to 



p(l) I pK^) I I ffV^) 



5(2) 



(1) I ^(2) 



-PT 



in accordance with (|3.4D. Including the ISR effects, we find that the expression 



(1,1,0) 



T2,max 



endpoint splits into two branches 



KX^LiMo,PT) = { 



F^^''''\m,,pt). 
fI^''^'\Mo,Pt) 



if Mo < Mo 
if Mo > Mo 



where 



(4.9) 
for the 

(4.10) 



4^'^'°)(Mo,Pt) 



/^(i,i,o)(pt) + 



/"(i,i,o)(Pr) + y) +M2 



(4.11) 



F^^'^'°)(Mo,pt) 



/^(i,i,o)(-Pt) + \l ( ^(1,1,0) (-pt) - ^ 



M^ 



and the py-dependent parameter /U(i^i^o)(pt) is defined as 

/^(l,l,0)bT) =/i(l,l,0) h/l+ 



PT 

2Mi 



Pt 



(4.13) 



Both branches correspond to extreme momentum configurations in which all three transverse 
vectors Pij^ , p^ and px are collinear. The difference is that the left branch F^'^'*^^ corre- 
sponds to the configuration i^'^ TT Pit ) TT PT, while the right branch f'^'^'^^ corresponds 
to (p^-p TT Pit ) l^J- Pt- Therefore, the two branches are simply related as 



F^^'^'°)(Mo,Pt) 



f1^'^'°)(Mo,-Pt). 



(4.14) 



It is easy to verify that in the absence of ISR, (i.e. for pT = 0) our general result ( |4.10| ) 



reduces to the previous formula ([4.31) . 



''This should not be confused with the term "balanced" used for the analytic Mt2 solutions discussed 
in [28,32]. 
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Figure 5: faj Dependence of the M^2^^l^ upper kinematic endpoint (solid lines) on the value of 
the test mass Mq, for Mi = 300 GeV, and Mq = 100 GeV, and for different values of the transverse 
momentum px of the ISR jet, starting from px = (green line), and increasing up to px = 5 TeV in 
increments of Apx — 100 GeV, from bottom to top. The uppermost line corresponds to the limiting 
case pT ^ oo. The horizontal (vertical) dotted line denotes the true value of the parent (child) mass. 
Solid (dashed) lines indicate true (false) endpoints. The red lines correspond to the function i^^^^^'"' 



defined in eq. (4.11 ), while the blue lines correspond to the function p^'^'*^^ dehned in eq. ( 4.1^ . (b) 



The value of the kink AO^^'i^o) defined in (4.21), as a function of the dimensionless ratios tt- 8,nd 



Mo 
Ah 



Our result (p^) for the A^li' max 

upper kinematic endpoint as a function of the test mass 
Mo is ihustrated in Fig. ||(a). We have chosen the same mass spectrum (Mq = 100 GeV and 
Ml = 300 GeV) as the one used in Ref. [31], so that our Fig. ^(a) can be directly compared 
to Fig. 9 of Ref. [31]. We consider a single ISR jet and show results for several different values 
of its transverse momentum pT, starting from pT = (the green solid line) and increasing 
the value of pT in increments of ApT = 100 GeV. The uppermost solid line corresponds 
to the limiting case Pt — > oo. The true value of the parent (child) mass is marked by the 
horizontal (vertical) dotted line. The red (blue) lines correspond to the function pj^'^'^^ 
^^a,i,o)^_ The solid port ions of those lines correspond to the true Myg^^^, endpoint, while 
the dashed segments are simply the extension of F^'^'^^ and F^'^'^^ into the "wrong" region 
for Mq, giving a false endpoint. 

Fig. ^(a) reveals that the two branches ( [4.11 ) and ( 4.12| ) always cross at the point 
(Mo,Mi), in agreement with eq. ( [4.71) . Interestingly, the sharpness of the resulting kink 
at Mo = Mo depends on the hardness of the ISR jet, as can be seen directly from ( 4.10 ). For 
small Pt, the kink is barely visible, and in the limit — > we obtain the old result ( [4.3| ) 
for the "balanced" momentum configuration, shown with the green solid line, which does not 
exhibit any kink. In the other extreme, at very large pT, we see a pronounced kink, which 
has a well-defined limit as pT oo. Our results in this regard are in agreement with the 
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findings of Ref. [31]. 

Tlie Af/^2^^ix kink exhibited in eq. ( |4l0| ) and in Fig. |(a) is our first, but not last, 
encounter with a kink feature in an M^g'^'^^ variable. Below we shall see that the Mt2 kinks 
are rather common phenomena, and we shall encounter at least two other kink types by the 
end of Sec. ^. Therefore, we find it convenient to quantify the sharpness of any such kink 
as follows. Consider a generic subsystem MljP^^'^^ variable whose endpoint Ml^^^^^r^iMc, pt) 
exhibits a kink: 



, , . . (4.15) 

4"'^''^(M„pt), ifM, >M,. 



The kink appears because -^r2'ma^a;(^c,PT) is not given by a single function, but has two 
separate branches. The first ("low") branch applies for Mc < Mc, and is given by some 

L 



function Fj^'^'^\Mc,PT)j while the second ("high") branch is valid for > M^, and is given 



by a different function, F^^'^''^\Mc,pt)- The function Mf^^^^axi-^^^ Pt) itself is continuous 
and the two branches coincide at Mc = M^- 

Ft^'''\Mc,PT) = F^^'^''\Mc,pt) , (4.16) 

but their derivatives do not match: 

(dFt!±\ , (dFf^ .4 

leading to the appearance of the kink. Let us define the left and right slope of the M!p2^^ax i-^c, Pt) 
function at Mc = Mc in terms of two angles and ' , correspondingly: 



tanei,"-'^-^^ ^{ '^^L y-^cj I ^ ^4_^g^ 




tan9i?'P'^^ = I , ^"^^^ I . (4.19) 



Mc=Mc 



Now we shall define the amount of kink as the angular difference AG^"''^''^^ between the two 
branches: 

tanG^ -tanG^ \ 

l + tanGi,"'^''=^tanGj^'^''V 

A large value of AG^"'^'"^-* implies that the relative angle between the low and high branches 
at the point of their junction Mc = Mc is also large, and in that sense the kink would be 
more pronounced and relatively easier to see. 

This definition can be immediately applied to the M^^^max kink that we just discussed. 
Substituting the formulas (4.11) and ( 4.121 ) for the two branches F^''^'^^ and F^'^'^^ into 
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the definitions (4.15,4.19) and subsequently into (4.20), we obtain an expression for the size 



AQi^'^fi) of the M?9^57kink: 



Mo (Mf - m2) pt ^ml+p^ 

Ml {Ml - M^f + 2M|f Ml {AMI + Pt) 



AG(^'1'°) = arctan ( '.,2,2^.1. ^^2. 1 • (4-21) 



The result ( 4.21 ) is illustrated numerically in Fig. ^(b). As can be seen from (|]2l|), Ae(^'^' 



L,0) 



depends on the two masses Mq and Mi, as well as the size of the ISR pT- However, since 
^0(1,1,0) jg dimensionless quantity, its dependence on those three parameters can be simply 
illustrated in terms of the dimensionless ratios ^ and This is why in Fig. |5|(b) we plot 
^0(1,1,0) degrees) as a function of ^ and 

Fig. |5|(b) confirms that the kink develops at large px, and is completely absent at = 0, 
a result which may have already been anticipated on the basis of Fig. |5|(a). For any given 
mass ratio the kink is largest for the hardest possible pT- In the limit pT ^ 00 we obtain 

r2 A/r2' 



hm Ae^^'°)= arctan , (4.22) 

PT^OO \^ 2M0M1 J ^ ' 



in agreement with the result obtained in [31]. From Fig. |5|(b) one can see that at sufficiently 
large pT, the AG^^'^'''^ contours become almost horizontal, i.e. the size of the kink Ae^^'^'") 
becomes very weakly dependent on pqp. A careful examination of the figure reveals that the 
asymptotic behavior at pT ^ 00 is in agreement with the analytical result ( [4.22| ). Notice 
that the maximum possible value of any kink of the type ( 4.15| ) is /S.Qmax^^ = 90°. According 



to Fig. 1(b) and eq. (^ ), in the case of AO^^'^'^) the absolute maximum can be obtained 



only in the — > 00 and Mq limit. The former condition will never be realized in 
a realistic experiment, while the latter condition makes the observation of the kink rather 
problematic, since the "low" branch Fl of the ^T'lmax^-^^iPT') function is too short to be 
observed experimentally. Therefore, under realistic circumstances, we would expect the size 
of the kink AG^^'^'*^-* to be only on the order of a few tens of degrees, which are the more 
typical values seen in Fig. |5|(b). 

According to Fig. ^(b), for a given fixed pt^ the sharpness of the AG'^^'^''') kink depends 
on the mass hierarchy of the particles Xi and Xq. When they are relatively degenerate, 
i.e. their iiictss ratio is large, the kink is relatively small. Conversely, when is much 

^ AT 

lighter than Xi, the kink is more pronounced. The optimum mass ratio which maximizes 
the kink for a given px, is rather weakly dependent on the pT, and for px ^ 00 eventually 



goes to zero, in agreement with eq. ( 4.22 ). However, for more reasonable values oi px as the 

Mo 
Ml 



ones shown on the left half of the plot, the optimal ratio ^ varies between 0.3 (at px ~ 0) 



to 0.1 (at Px ~ 5Mi). In this sense, the value of ^ = | which was chosen for the illustration 
in Fig. |5|(a) (as well as Fig. 9 in Ref. [31]) is rather typical. 

In conclusion of this subsection, it is worth summarizing the main points from it. The 
good news is that the AG'-^'^''^^ kink in principle offers a second, independent piece of in- 
formation about the masses of the particles Xq and Xi. When taken together with the 
^T2^max eudpoiut measurement ( |4.5| ), it will allow us to determine both masses Mq and Mi, 
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in a completely model-independent way. Our analytical results regarding the AO^^'^'''^ kink 
complement the study of Ref. [31], where this kink was first discovered. However, on the 
down side, we should mention that much of our discussion regarding the AQ^^'^'^^ kink may 
be of limited practical interest, for several reasons. First, as seen in Fig. ^, the kink becomes 
visible only for sufficiently large values of the pT- Since the ISR px spectrum is falling rather 
steeply, one would need to collect relatively large amounts of data, in order to guarantee the 
presence of events with sufficiently hard ISR jets. Even then, the collected events may not 
contain the momentum configuration required to give the maximum value of M^^^'°\ An 
alternative approach to make use of the kink structure would be to measure the endpoint 
function -^^^'^^^^(-^OiPt) for several different px ranges, and then fit it to the analytical 
formula ( |4.10D . Whether and how well this can work in practice, remains to be seen, but the 
results of [31] from a toy exercise in the absence of any backgrounds and detector resolution 
effects do not appear very encouraging. Nevertheless, while the kink structure A©^^'^''^) may 
be difficult to observe, the measurement ( |4.5| ) of the endpoint -^r2maa;(-^o ~ 0,pT = 0) 



should be relatively straightforward. In Sees. |5.1| and 5.3 we shall see that the additional 
Mt2 information from events with n = 2 decay chains will eventually allow us to determine 
all the unknown masses. 

4.2 The subsystem variable M^g^'^^ 

The subsystem variable M^g'^'^^ is illustrated in Fig. §(b), where we use the subchain within 

(2 2 1) 

the smaller rectangle on the left. is a genuine subchain variable in the sense that we 

only use the SM decay products X2, and ignore any remaining objects arising from the two 
xi's. In the absence of ISR {pT = 0) one can adapt the results from [32] and show that the 

(2 2 1) 

formula for the Mrp2 ' endpoint is 



M^'i^LiMuPT = 0) = ^(2,2,1) + ^4,2,1)+^? , (4.23) 



where the parameter ^(2,2.1) was defined in eq. (|4 



_M2( Ml\_ M^-Ml 



Almost all of our discussion from the previous Section 4.1 can be directly applied here as 
well. For example, in order to measure the parameter /i(2,2,i)) we only need to extract the 
endpoint of a single distribution, for a single fixed value of the test mass Mi . As before, we 
choose to use Mi = 0. The resulting endpoint measurement 

M^SLiMi =0,pT = 0) = 2/.(2,2,i) = -^M^ (4-25) 



provides the required measurement of the parameter /U(2,2,i) appearing in eq. ( 4.23| ), as well 



as one constraint on the masses Mi and M2 involved in the problem. More importantly, the 



new constraint ( 4.25| ) is independent of the relation ( [4.5D found previously in Sec. 



(221) (221) 
The new variable ' will also exhibit a kink in the plot of its endpoint Mrj.2max 

as a function of the test mass Mi. This is the same type of kink as the one discussed in 
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the previous subsection, therefore all of our previous results would apply here as well. In 
particular, the analytical expression for the kink is given by 



( ^^2,2,1)^^^^^^^ ^ if Ml < Ml 

^Fg'^'^\m,PT), if Ml > Ml 



where 



/^(2,2,l) (Pt) + \l ( /^(2,2,1) (Pt) + y) ^1 



1 

2 ^2 

Pt 



/^(2,2,i)(-Pt) + \ (/^(2,2,i)(-Pr) 



PT 



Pt 



(4.26) 



(4.27) 
(4.28) 



and the pr-dependent parameter /^(2,2,i)(pt) is defined in analogy to (|4.13[) 



^(2,2,1) (pt) = /^(2,2,1) \/l + 



PT 
2Mi 



PT 

2Mo 



(4.29) 



The size of the new kink Ae(2'2.i) 

can be easily read off from eq. (|4.21[) , where one should 
make the obvious replacements Mq — > Mi and Mi M2. 

We can now generalize the two examples discussed so far (M^g^''^^ and M-^2'^'^^) to the 
case of an arbitrary grandparent index n, with p = n and c = n — 1. We get 



M 



{n,n,n—l) 



T2,max 



(Mn-UPT) 



p{n,n,n—l) 



F<f'"'"-^)(M„_i,pT) , if M„_i > M„_i , 



(M„_i,pt), ifM„_i<M„_i, 



(4.30) 



where 

jp{n,n,n—l) 



{Mn-l,PT) = { 



n,n,n—l) (pt) + \ (PT) + f) +MLl 



(n,?i,n— 1) / X 
^ '{Mn-1,PT) 



f^{n,n,n-l)i-PT) + \ ( /^(n,n,n-l) ("Pt) 



1 

2 2 ^ ^ 

Pt 



(4.31) 



4 

(4.32) 



and the pT-dependent parameter /i(n,n,n-i)(PT) is simply the generalization of eqs. ( 4.13 ) and 



(pt) M(n,n,?i— 1) 



'1 + 



PT 



2M„/ 2M„ 



(4.33) 



For n = 1 or n = 2, the general formula ( [4. 301 ) reproduces our previous results ( 4.10| ) and 
(4.26), correspondingly. 
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('2 2 0) 

4.3 The subsystem variable Mj,^ 

The variable M/j^g'^'''^ is illustrated in Fig. Q(b), where we use the whole chain within the 
larger rectangle. As long as we ignore the effects of any ISR, we have a balanced^ momentum 
configuration and the analytical results from Ref. [32] would apply. In particular, the endpoint 
M^X^LiMo,PT = 0) is given by [32] 

(2 2 0).~ . f Ff''°)(Mo,PT = 0), if Mo < Mo, 

<2±(Mo,PT = 0)= ; ^ (4.34) 

[ F^''2'°) {Mo,PT = 0) , if Mo > Mo , 

where 



F^^^'''\Mo,Pt = 0) = ^(2,2,0) + ^4,2,0) + ^o> (4.35) 



F^'^'^\Mo,Pt = 0) = /i(2,2,l) + M(2,l,0) + \/ {^^{2,2,l) - /^(2,1,0))^ + M^ , (4-36) 

and the various parameters H{n,p,c) a-re defined in ([4.1|). Notice that these expressions are 
valid only for px = ^- We have also derived the corresponding generalized expression for 
M^2'^aa,(-Mo,pT) for arbitrary values of pT, which we list in Appendix^ 

The most striking feature of the endpoint function ( [4.34 ) is that it will also exhibit a 



kink AG^-^'^'*^) at the true value of the test mass Mo = Mo. However, as emphasized in [31], 
the physical origin of this kink is different from the kinks AB^^'^''') and AO^^'^'^^ which we 



encountered previously in Sees. 4.1 and 4.2. This is easy to understand - in Sees. and 4.2 



we saw that the kinks Ae(^'^'°) and AO^^'S.i) 

arise due to ISR effects, while eq. ( |4.34D holds 



in the absence of any ISR. The explanation for the A©^'^'^''^) kink has actually already been 
provided in [32]. In essence, one can treat the SM decay products xi and X2 in each chain as 
a composite particle of variable mass, and the two branches pj^''^'^^ and F^'"^'^^ correspond 
to the two extreme values for the mass of this composite particle. 



In spite of its different origin, the kink in the function (4.34) shares many of the same 
properties. Let us use a specific example as an illustration. Consider a popular example from 
supersymmetry, such as gluino pair-production, followed by sequential two-body decays to 
squarks and the lightest neutralinos. This is precisely a cascade of the type n = 2, in which 
X2 is the gluino ^, Xi is a squark q, and Xq is the lightest neutralino Xi- Let us choose the 
superpartner masses according to the SPSla mass spectrum, which was also used in Ref. [32]: 

M2 = 613 GeV, Ml = 525 GeV, Mo = 99 GeV. (4.37) 

The resulting function M^2'^ix(Mo,PT = 0) is plotted in Fig. |6|(a) with the upper set of lines 
(compare to Fig. 12(b) in Ref. [32]). 

{'220) ^ 

There are several noteworthy features of M^g ^aa;(Mo,PT = 0) which are evident from 
Fig. ^(a). First, when the test mass Mq is equal to the true child mass Mo, the Mt2 endpoint 
yields the true parent mass, in this case M2: 

M^f^LiMo = Mo,PT = 0) = M2 . (4.38) 



In the sense of eq. (4.8). See the discussion following eq. 
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Figure 6: Dependence of the Mrp^ ^ and Mj^^ upper kinematic endpoints on tie value of the 
test mass Mq, for (a) the SPSla parameter point in MSUGRA: Ah = 613 GeV, Mi = 525 GeV, and 
Mo = 99 GeV; or (b) a spht spectrum M2 = 2000 GeV, Mi = 200 GeV, and Mq = 100 GeV. The 
horizontal (vertical ) dotted lines denote the true value of the parent (child ) mass for each case. Solid 
(dashed) lines indicate true (false) endpoints, while red (blue) lines correspond to F^'P''^^ ^p^^P^'^) ^ 
branches. 



This property of Mj'2 is true by design, and is confirmed by the dotted hnes in Fig. |6|(a). 

I 1 (2 2 0) ~ 

Second, as seen from eq. ([4.34|) , Mj,2 max(-^o,PT = 0) is not given by a single function, but 

(2 2 0) ~ 

has two separate branches. The first ("low") branch Fj^' ' ' apphes for Mq < Mq, and is 

_ (2 2 0) 

shown in Fig. y(a) with red hnes. The second ("high") branch ' is valid for Mq > Mq 
and is shown in blue in Fig. ^(a). While the two branches coincide at Mq = Mq: 

Ff '''°)(Mo,pT = 0) = 4''''°)(Mo,PT = 0) , (4.39) 
their derivatives do not match: 

'2.0) \ /;5p(2.2,o)\ 

(4.40) 

Mo=Mo 

leading to a kink Ae^^.s.o) in the function M^^^^I^{Mq,pt = 0) [29-32]. As before, let us 
try to investigate quantitatively the size of this kink. Applying the general definition ( ^.2[1| ), 
we obtain 




Ae(2'2'0) = arctan ( - y)(^ - ^)Vf ) , (4.41) 
\{y + z){l+yz)+4yzl ^ ^ 



where we have defined the squared mass ratios 

Ml ^ Ml 

ir Ml 



y = ^, z = ^. (4.42) 

^ Ml M? ^ ^ 



The result (4.41) is plotted in Fig. ^(a) as a function of the mass ratios and ^/z. Fig. |^(a) 



demonstrates that as both y and z become small, the kink A0*^2>2,o) g^^g ^lore pronounced. 
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Figure 7: The amount of kink: (a) AO*^^'^^*^) and (b) AB*^^'^^"-' in degrees, as a function of the mass ra- 
tios and \fz. The white dot and the white asterisk denote the locations in this {^/y, \fz) parameter 
space of the two sample spectra (4.37 ) and {4.5E ) used for Figs, ^(a) and ^(b), correspondingly. 



Fig. I^a) also shows that the kink AG^'^'^''^^ is a symmetric function of y and z, as can also be 
seen directly from eq. (i-41). Therefore, the kink AB^^'^''^^ will be best observable in those 
cases where y and z are both small, and in addition, the mass spectrum happens to obey the 
relation'' y = z, i.e. Mi = \JMqM2. Unfortunately, the SPSla study point is rather far from 
this category - the spectrum ( |4.37] ) corresponds to the values y/y = 0.856 and ^/z = 0.189, 
which are indicated in Fig. ^(a) by a white dot. This conclusion is also supported by Fig. |6|(a), 
which shows a rather mild kink in the SPSla case. 

We shall be rather ambivalent in our attitude toward the AO'^^'^''') kink as well. While the 
interpretation of the kink is straightforward, its observation in the actual experiment is again 
an open issue. On the one hand, the experimental precision would depend on the particular 
signature, i.e. the type of the SM particles xi and X2- If those are leptons, their 4-momenta 

(k) (k) 

PI and P2 will be measured relatively well and the kink might be observable. However, 
when xi and X2 are jets, the experimental resolution may not be sufficient. Secondly, as seen 
in Fig. ^(a), the kink itself may not be very pronounced, and its observability will in fact 
depend on the particular mass spectrum. 

The main lesson from the above discussion is that while the existence of the kink is 
without a doubt, its actual observation is by no means guaranteed. Therefore, our main 



mass measurement method, described later in Sec. 5.1, will not use any information related 



to the kink. In fact in Sec. |5.1j we shall show that one can completely reconstruct the mass 
spectrum of the new particles, using just measurements of Mt2 endpoints, each done at a 
single fixed value of the corresponding test mass. It is worth noting that, in general, an 
endpoint in a spectrum is a sharper feature than a kink of the type ( 4.20| ). Therefore, we 



^Notice that for this special value of Mi — \/MoM2, the upper endpoint of the invariant mass distribution 
Mxix2 is the same as in the case when the intermediate particle Xi is off-shell, i.e. when Mi > M2. Then we 



find that the M 



(2,2,0) 
T2,max 



cases. For details, see Appendix ^. 



formulas and corresponding kink structures are identical in the on-shell and off-shell 
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would expect that the experimental precision on the extracted endpoints will be much better 
than the corresponding precision on the kink location. The kink will also not play any role 
in our hybrid method, described in Sec. Only for the method described in Sec. |5.2| , we 
shall try to make use of the kink information. 



Let us now return to our original discussion of the ' ' endpoint (4.34). Following our 



previous approach from Sees. 4.1 and 4.2, we would choose a fixed value of the test mass M( 







and measure the corresponding Mt2 endpoint. However, the presence of two branches ( 4.35| ) 



and ( [4.36 ) leads to a slight complication: for a randomly chosen value of Mq, we will not 



know whether we should use ( 4.351) or ( 4.36| ) when interpreting the endpoint measurement. 



This requires us to make very special choices for the fixed value of Mq, which would remove 
this ambiguity. It is easy to see that by choosing Mq = 0, we can ensure that the endpoint 
is always described by the "low" branch ( [4.35| ), and the Mj^^max measurement can then be 
uniquely interpreted as 

<2il (Mo = 0, = 0) = 2 ^(2,2,0) = • (4.43) 

However, we could also design a special choice of Mq, which would select the "high" branch 



(4.35) and again uniquely remove the branch ambiguity. For this purpose, we must choose a 
value for the test mass Mq which is sufficiently large, in order to safely guarantee that it is 
well beyond the true mass Mq. Since the true mass Mq can never exceed the beam energy 
Ef,, one obvious safe and rather conservative choice for Mq could be Mq = Ef^, in which case 
from ( 4.34| ) we get 



M&mi(Mo = Eb,PT = 0) = /i(2,2,l) + /X(2,l,0) + V (^(2,2,1) " /^(2,1,0))' + E^ (4.44) 



M2/Mf M|\ Ml ( Ml _ Ml 
2 MfJ^y 4 \Ml Ml I ^ 



Notice that the high branch function ' in eq. (4.36) is rather unique in one very impor- 



tant aspect: it depends not just on one, but on two mass parameters, namely the combinations 

f2 2 0) 

/^(2,2,i) + /"(2,i.o) s-'^d /^(2,2,i) ~ /^(2,i,o)- Contrast, the "low" branch ' , as well as the 
previously discussed endpoint functions Af^g'^ax ^^(MojPt = 0), each contained a single /i 
parameter. As a result, in those cases we did not benefit from any extra measurements for 
different values of the test mass Mq ~ had we done that, we would have been measuring the 

(2 2 0) 

same ^ parameter over and over again. However, the situation with F^ ' is different, and 
here we will benefit from an additional measurement for a different value of Mq. For example, 
let us choose Mq = E'^, with E"^ > Ef,, which will still keep us on the high branch. We obtain 
another constraint 



KXrnU^O = EipT = 0) = /i(2,2,l) + A^(2,l,0) + V (/^(2,2,1) " M(2,l,0))' + ^fe' (4-45) 



M2 (Ml Mi\ MUMl M2 



1 2 
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It is easy to check that the constraints (4.45- 05| ) are all independent, thus providing three 
independent equations^ for the three unknown masses Mq, Mi and M2. These three equations 
( |4.43 - 44^ ) can be solved rather easily^, and one obtains the proper solution for the masses 
Mo, Ml and M2, up to a two-fold ambiguity: 



M2 ^ M2 



Ml 



Mo 
Ml 



M2 , Mo ^ Mo 



(4.46) 



which is nothing but the interchange y <-> z at a fixed M2. The ambiguity arises because the 
expression ( [4.34 ) for the endpoint M,|?2^^^ (and consequently, the set of constraints ( [4.43| - 
4.45| )) is invariant under the transformation (4.46). Because of this ambiguity, in addition to 



the original SPSla input values (4.37) for the mass spectrum, we obtain a second solution 



M2 = 613 GeV, Ml = 115.6 GeV, Mq = 99 GeV. 



(4.47) 



This second solution was missed in the analysis of Ref. [32]. It is easy to check that the 
alternative mass spectrum ([4.47] ) gives an identical j^^^{Mq,pt = 0) distribution as 

the one shown in Fig. |6|(a), so that it is impossible to rule it out on the basis of Mj^^^^^ 
measurements alone. 

The previous discussion reveals an important and somewhat overlooked benefit from 
the existence of the kink - one can make not one, not two, but three independent endpoint 
measurements from a sing le m!^/'''^ distribution! In fact, we shall argue that the three 
measurements ( [4.43 -4.45) are much more robust than the kink measurement ( [1.2C1[ ). For 



(2 2 0) 

example, when the child mass is relatively small, the lower branch i*^ ' is relatively short 
and the kink will be difficult to see, even under ideal experimental conditions. An extreme 
example of this sort is presented in Section ^, where we discuss top quark events, in which the 
child (neutrino) mass Mo is practically zero and the kink cannot be seen at all. However, even 
under those circumstances, the endpoint measurements ([4.43 - 05| ) are still available. More 



importantly, the constraints (4.43-4.45) are independent of the previously found relations 
( [4.5] ) and ( [4.25 ), so that the latter can be used to resolve the two- fold ambiguity ( [4.46| ). 

Before we move on to a discussion of the last remaining subsystem Mt2 quantity in the 



M^g'^'*^^ variable yields three independent endpoint measurements (4.43-4.45), and possibly a 



next Sec. 4.4, let us recap our main result derived in this subsection. We showed that the 



kink measurement ( 4.2C| ). The M^g'^''''* endpoint measurements by themselves are sufficient 



to determine all three masses Mq, Mi and M2, up to the two- fold ambiguity ( [4.46 ). This 
represents a pure MT2-based mass measurement method, which does not use any any kink 
or invariant mass information. 



*In practice, instead of relying on individual endpoint measurements for three different values of Mq, one 
may prefer to use the experimental information for the whole function ^^^^^^^.(A^OjPt ~ 0) and simply fit to 
it the analytical expression (4.34) for the three floating parameters Mo, Mi and M2, as was done in Ref. [32]. 
As we shall see shortly, this method does not lead to any new information, and may only improve the statistical 
error on the mass determination. Therefore, to keep our discussion as simple as possible, we prefer to talk 
about the three individual measurements (1.43 - 4^4^ ) as opposed to fitting the whole distribution (4.34). 

^The general solution for M2, Mi and Afo in terms of the measured endpoints ( 4.43| -4.45) is rather messy 
and not very illuminating, therefore we do not list it here. 
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4.4 The subsystem variable m!^?^^'^^ 

The variable is illustrated in Fig. ^(b), where we use the subchain within the smaller 

rectangle on the right. This is another genuine subsystem quantity, since we only use the 
SM decay products xi and ignore the upstream objects X2- However, the upstream objects 
X2 are important in the sense that they have some non-zero transverse momentum, and as 
a result, the sum of the transverse momenta Pqj. of the children Xq is not balanced by the 
sum of the transverse momenta of the SM objects xi used in the Mt2 calculation: 



5(1) 



5(2) 



(1) I ^(2) 



P2T 



-(2) 

- 



PT^O . 



(4.48) 



Notice that even in the absence of any ISR px, this is still an unbalanced configuration, due 
to the transverse momenta P2T upstream objects X2- Therefore, we cannot use 

the existing analytical results on Mt2, since previous studies always assumed that the right- 
hand side of eq. ( 4.481) is exactly zero, due to the lack of any particles upstream. We therefore 
need to generalize the previous treatments of Mt2 and obtain the corresponding endpoint 
formulas for our new subsystem m!^^^'^^ variable. The general mathematical properties of 
subsystem M!j?2^'^^ variables will be presented in a forthcoming publication [49]. Here we 
shall only use the results relevant for our example (n = 2). In particular, in the absence of 



any intrinsic ISR (i.e., for pT 
given by 



(2,1,0) 



T2,max 



(Mo,PT = 0) 



0), we find that the endpoint of the 

7(2,1,0) 



(2,1,0) 



distribution is 



R 



'{Mo,pT 



0), 
0), 



if Mo< Mo, 
if Mq> Mo, 



where 



p(2,l,0) 



{Mo,PT = 0) 



/^(2,2,0) - M(2,2,l) + 



/^(2,2,0)+^0 



2 

^(2,2,1) 



(4.49) 



(4.50) 



7^(2,1,0) 
R 



{Mo,PT = 0) 



M(2,l,0) + V (^(2,2,1) - /^(2,1,0))' 



+ M2 



^(2,2,1) 



(4.51) 



and the various parameters fJ'(^n,p,c) are defined in (4.1 



general pT (i-e., arbitrary intrinsic ISR) are listed in Appendix]^. 

From eq. ( ^.49| ) we see that, once again, the endpoint function m!^2^^1x(^o^Pt 
would exhibit a kink AG^^'^'^^^ at the true value of the test mass Mq = Mo- 



The corresponding expressions for 

0) 



(2,1,0)' 



dMo 



(2,1,0)' 



R 



dMn 



(4.52) 



A/o=Afo 

The existence of this kink should come as no surprise 



' Mo=Mq 

Ref. [31] showed (in the px 



00 



limit) that any type of upstream momentum will generate a kink in an otherwise smooth 
Mt2. max function. As before, the value of the M'j'2 endpoint -/Vf^2 ^Ix 
reveals the true mass of the parent: 



at the kink location 



(2,1,0) 
T2,max 



(Mo = Mo,PT = 0) = Mi 



(4.53) 
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At the same time, the physical origin of this kink is different from either of the two kink types 
(/X0(i.i:O) /\0(2,2,o)-j (Jigcussed earlier. Clearly, the new kink is different from AQ^'^''^'^\ 
which was due to the varying invariant mass of the {xi,X2} system. Here we are using a 
single SM particle xi whose mass is constant. Furthermore, the new kink Ae(2.i'0) 

cannot 

be due to any ISR like in the case of A©^^'^''^), since eq. ( [4.49 ) does not account for any ISR 



effects. The real reason for this new is a third one, namely, the kinematical 

restrictions placed by the decays of the upstream particles (in this case, the grandparents 
X2). 

We now proceed to investigate the new kink AO^^'^''^^ quantitatively. Using the same 



example of gluino pair-production for the SPSla mass spectrum (|4.37D , we plot the function 



( |4.49 ) in Fig. ^(a). Comparing the lower and the upper set of lines in the figure, we notice 



that the m!^2^^1x ^^'^ ^T2^max variables share several common characteristics. They both 
exhibit a kink at the true location of the child mass Mq = Mq, while their values at that 
point reveal the true parent mass in each case: Mi for Mf^^^^lx ^^'^ -^2 for M^^'^^^^. Using 
the definition ( ^) , we find that the size of the Ae(2.i'0) kink is given by 



Ae(^-^-°) = arctan ( ^ , ) , (4-54) 



where the parameters y and z were aheady defined in ( |4.42D . The kink AG^^.i-O) plotted 



in Fig. |^(b) as a function of ^Jy and y/z. We notice that the kink structure becomes more 
pronounced for relatively small y and z. Comparing Figs. |^(a) and|^(b), we see that for any 
given set of values for y and z, the AG^^'^'^^^ kink discussed here is more pronounced than the 
Ae(2,2,o) Yxv^ from the previous subsection^'^ . The difference is particularly noticeable in the 



region of ~ and \fz ~ 0.2. The SPSla mass spectrum ( 4.37 ) in our previous example 
was rather far away from this region, as indicated by the white dots in Fig. ^. Now let us 
choose a different mass spectrum, which is closer to the region where the difference between 
the two kinks becomes more noticeable, for example 

M2 = 2000 GeV, Ml = 200 GeV, Mq = 100 GeV, (4.55) 

corresponding to the point marked with the white asterisk in Figs. |^(a) and |^b). The 
resulting endpoint functions M^^max ^T2max plotted in Fig. ^(b). Indeed we see 
that with this new spectrum the kink in the ^'^2max function is much more noticeable than 

the kink in the M/^g'^Qj, function. Therefore, our first conclusion regarding the m}y2'^^ 

variable is that its kink is in general sharper and appears to be more promising than the 

{1 2 0) I — 
previously discussed kink in the Mj.^ ' ' variable from Sec. 4.3. 



Following our previous strategy, we shall not dwell too long on the kink, but instead 
we shall discuss the available endpoint measurements for various values of Mq. Again, the 
presence of two branches in eq. ( 4.49| ) can be used to our advantage. As in Sec. 4.3, we first 



choose a test mass value Mq = 0, which would "select" the low branch ( [4. 501 ) and result in 



This statement can also be verified using the analytical formulas (4.41) and (tt.S- 
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an endpoint measurement 



0,PT = 0) = 2 a/ M(2,2,0) (^(2,2,0) - /^(2,2,1)) • 



(4.56) 



Just as before, we could also choose a rather large value for Mq 
high branch ( [4.51| ) and result in the measurement 



Eb, which would select the 



Mr 



(2,1,0) 



T2,max 



/^(2,l,0) + Y (/^(2,2,1) - /^(2,1,0)) + 



^(2,2,1) 



(4.57) 



A third choice, Mq = Ei, with E', > E^, would yield yet another endpoint measurement 



(2,1,0) 



T2,max 



/^(2,l,0) + V (^(2,2,1) - /^(2,1,0)) + E't 



M^2,2,l) 



(4.58) 

Again we obtained three equations (4.56- 4?5^ ) for the three unknown /i-parameters /i(2,2,o); 
/^(2,2,i) M(2,i,o)5 equivalently, for the three unknown masses Mq, Mi and M2. These 
equations are all independent and can be easily solved, giving a total of four solutions. 
However, three of the solutions are always unphysical, so that we end up with a single unique 
solution. This represents an important advantage of the Af^2^ix variable in comparison 



(2 2 0) 

with the Mrp^'' variable discussed in Sec. 



There we found that M, 



(2,2,0) 



T2,max 



always gives 

rise two a two-fold ambiguity in the mass spectrum, while now we see that Mf^^^^lx does 
not suffer from this problem and already by itself allows for a complete and unambiguous 
determination of the mass spectrum. 



5. Mr2-based mass measurement methods 

In this section we use the analytical results derived in the previous section to propose three 
different strategies for determining the masses in n < 2 decay chains. We shall illustrate 
each of our methods with a specific example, for which we choose to consider the dilepton 
samples from W~^W~ and tt events. The former is an example of the n = 1 decay chain 
exhibited in Fig. §(a), while the latter is an example of the n = 2 decay chain in Fig. §(b). 
Most importantly, these samples already exist in the Tevatron data and will also be among 
the first to be studied at the LHC. Correspondingly, throughout this section we shall use the 
following mass spectrum 

M2 = mt = 173 GeV, 

Ml = mw = 80 GeV, (5.1) 
Mo = =0 GeV. 

Before we begin, let us review the four different iW^g'^'^^ variables which are in principle 
available in that case. Each one of them is plotted in Fig. R for five different values of the 
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Figure 8: Unit-normalized distributions of M, 



(n,p,c) 
T2 



variables in dilepton events from (a) W^W^ 
pair production and (b-d) tt pair production. Each panel shows results for five different values (0, 
100, 200, 300 and 400 GeV) of the corresponding test mass. The methods of Sees. 5.1 and only- 
make use of the Mt2 endpoint at zero test mass, M^^'^^^^{Mc = 0), which is indicated by the vertical 

red arrow. In panel (c), the two dotted line M^j'^'"' distributions correspond to the correct and the 
wrong pairing of the two b-jets with the leptons, while the solid line distribution is the average of 
these two. 



corresponding test mass (0, 100, 200, 300 and 400 GeV). In Fig. |(a) we show the All^f'^' 
variable from W^W~ pair production events, while in Figs. §(b-d) we correspondingly show 
the M^2'^'^\ m!^^'^'^^ and M^jp^^'^^ variables from ti events. We used PYTHIA for event 
generation and did not impose any selection cuts, since they will not affect the location of 
the Mt2 endpoint^^. The plots are made for the Tevatron (a pp collider with a 2 TeV center- 
of-mass energy), where the relevant data is already available. The corresponding analysis 
for the LHC is very similar. All of our plots in this Section have the full ISR effects. As 
discussed in Section ^ the presence of ISR with nonzero px will increase the nominal M^g'^'*^^ 

^^The cuts would have an impact on the overall acceptance and efficiency. This eflFect is not relevant here, 
since we are showing unit-normalized distributions. The cuts may also distort the shape of each distribution, 
but should preserve the location of the upper endpoint. 
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endpoints: 

> M^^^{M^,&) , (5.2) 

where the equahty is obtained only when Mc = Mc- ISR will therefore introduce some 
systematic error when one is trying to measure ^j^^axi^f:^ '^^^ ^^^^ ^^^^ error depends 
on the ISR pt spectrum, which in turn depends on the type of collider (Tevatron or LHC). 
At the Tevatron, this will not be such a serious issue, as evidenced from Fig. ^, where the 
observed endpoints in the presence of ISR match pretty well with their expected values for 
the Pt = case. On the other hand, at the LHC this may become a problem, which can 
be handled in one of two ways. First, depending on the particular signature, one may be 
able to select a sample with « (at a certain cost in statistics), by imposing a suitably 
designed jet veto to remove jets from ISR. Alternatively, one can use the full event sample 
(which would include ISR jets), and make use of our general formulas in Appendix |A|, which 
contain the explicit pT dependence of M^2^^1x- 

In the previous Section ^ we derived that in the case of n < 2 cascades, there are 8 
different Mt2 endpoint measurements: one for M^2^^lx (^^^ 14-1|) , one for 

M^2' ml (see eq. ( pSD and Sec. |4|), three for M^g'm^l (see eqs. (|03[^l5| ) and Sec. |]D, and 



three for m|?2^^^^ (see eqs. ( [4.56| - ^l58[ ) and Sec. iA). Given that we are trying to determine 



only three masses Mq, Mi and M2, it is clear that these 8 measurements should be sufficient 
to completely determine the spectrum. We also see that our previous count (|3.8D has actually 
greatly underestimated the power of M'r2, and the number of available measurements is in 



fact much larger than the number of M^f^^^^ variables. Indeed, as shown in Sections O 



and 4.4, there are cases where we might be able to obtain more than one mass constraint 



from a given M^f^^^^ variable. Of course, the 8 measurements cannot all be independent 
among themselves, as they only depend on three parameters. Our three methods below will 
be distinguished based on which subset of these measurements we are using. 

5.1 Pure Mt2 endpoint method 

With this method, we use endpoint measurements Enpc at a single fixed value of the 
test mass, which for convenience we take to be Mc = 0: 

Enpc ^ M^^^^l {Mc = 0,PT = 0). (5.3) 

The corresponding formulas interpreting those measurements in terms of the physical masses 
Mq, Ml and M2 were derived in Sec. 

^110 - <2'±(o,o) = = M2Vy{i - z) , (5.4) 

^221 ^ M^'^liO, 0) = ^^^^ = M2{l-y), (5.5) 
E220 - M^^iijO, 0) = = M2{1- yz) , (5.6) 



£'210 = 



M2 

Sil(0, 0) = {Ml - M^mf - M^) = M2Vy{l-z){l-yz) . (5.7) 
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Using the mass spectrum ( |5.1| ), the predicted locations of these four Mt2 endpoints are 



-^110 


= 80 GeV , 


-E221 


= 136 GeV 


-£■220 


= 173 GeV 


-£^210 


= 80 GeV , 



(5.8) 
(5.9) 
(5.10) 
(5.11) 



which are marked with the vertical red arrows in Fig. |^. Given that we have four measure- 
ments ( |5.4f|5.7| ) for only three parameters Mq, Mi and M2, one should be able to uniquely 
determine all three of the unknown parameters. Naively, it seems that using just three of the 
measurements ([5. 4-5. 7) should be sufficient for this purpose, and furthermore, that any three 
of the measurements (5.4-^?7) will do the job. However, one should exercise caution, since 
not all four measurements (^-5.7) are independent. It is easy to check that £'221, -E220 and 
-E210 obey the following relation 



-^210 — -^220 (£-220 — -E'221) 



(5.12) 



This means that in order to be able to solve for the masses from eqs. (5.4-x^), we must always 
make use of the Ehq measurement in eq. (5.4), and then we have the freedom to choose any 
two out of the remaining three measurements ( |5.5| - p.7] ). For example, using the set of three 
measurements {-Eno, -£'2215 -£'220} (i-e- eqs. (^.4H5.6|) ), the masses are uniquely determined as 



Mn 



Ml 



Mo 



-E'lio {-E'221 (-6^220 — -£-221) [-E'220 (-£-220 — -£-221) — -E-iio] }■ 



-E'lio -E'221 {E; 



220 



^110 ~ (^220 
- -E'221) 



-£■221)^ 



-^110 



(i?220 — -E'221)^ 

EiiQ — {E220 — E221Y 



(5.13) 
(5.14) 
(5.15) 



Similarly, one can solve for Mq, Mi and M2 using the set of measurements {-Eno, -£"220, -E'210}) 
or alternatively, the set of measurements {-Eno, -©221, -E'2io}- In each case, the remaining 
fourth unused measurement provides a useful consistency check on the mass determination. 



5.2 M7^2 endpoint shapes and kinks 



The method proposed in Sec. 5.1 uses the measured endpoints from several different Mrp^ 



{n,p,c) 



variables. Now we discuss an alternative method which makes use of a single M 



(n,p,c) 
T2 



variable. 



Let us begin with the simplest case of n = 1 as shown in Fig. ^a). In that case, we 



have only one Mt2 variable at our disposal, namely M 

r(l:l.O) " 



(1,1,0) 
T2 



Its properties were discussed 

in Sec. where we showed that its endpoint Mj^g^^jj, can allow the determination of both 



masses Mq and Mi, at least as a matter of principle. Indeed, the endpoint measurement 
( |5.4| ) at zero test mass provides one relation among Mq and Mi. The key observation in 
Sec. [4.1| (which was first done in [31]) was that with the inclusion of ISR effects, the endpoint 
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function M^^^^^^{Mq, pt) exhibits a kink at Mq = Mq, which can then be used to determine 
both masses Mq and Mi. The method can be readily apphed to the existing dilepton event 
sample from W~^W~ pair production, which will allow an independent measurement of the 
W mass mw and the neutrino mass rrii^. While the precision of this measurement will not 
be competitive with existing W and neutrino mass determinations, it is nevertheless useful 
to test the viability of this approach with real data. 

Now let us discuss the more complicated case of n = 2, which in our example corresponds 



to tt pair production with both tops decaying leptonically. As discussed in Sees. 4.2, 4.3 and 



4.4| , here we have a choice of three different Mt2 variables: Mj^^ ' ■> ' ■> and Mrp^ 



(2,2,1) ,^(2,2,0) . ,712,1,0) 



Because of the larger tt cross-section, we expect that the statistical precision on each one 
of those three variables will be better than the M^g'^''^'' variable of the n = 1 case. As 



f 2 2 0) ("2 1 O") 

shown in Sees. [4.3| and |4.4|, each of the two variables ' and iVf^g ' exhibits a kink 



in its endpoint Ml^^^^ax '^^en considered as a function of the test mass Mq, even when the 
transverse momentum of the intrinsic ISR in the event is zero, pT = 0. Then, which of these 

f2 2 0) ~ 

two variables is better suited for a mass determination? The case of Mj.2 maxi^o^PT = 0) was 
already discussed in [29,31,32]. Here we would like to propose the alternative measurement 
of M!j^2^^lx(^o^PT ~ What is more, we would like to emphasize that our function 
^T2mi(^o,PT = 0) offers several unique advantages over the previously considered case of 

1. The subsystem variable M^g'^'*^^ does not suffer from the combinatorics problem which 
is present for M!jp2^'^\ Indeed, when constructing the M^2^'^^ distribution, one has 
to decide how to pair up the 6-jets with the two leptons. Because it is difficult to 
distinguish between a b and a b, there is a two-fold ambiguity which is quite difficult 
to resolve by other means. In contrast, our subsystem variable M!jp2^'^^ does not make 
direct use of the 6-jets, and is therefore free of such combinatorics issues. 



2. As we already saw in Sec. 4.2, even under perfect experimental conditions, the fit to the 



(2 2 0) 

Mrp2 ^^2, endpoint results in two separate solutions for the mass spectrum: one solution 



(see ( ^.37 )) is given by the true values of the input masses, while the second solution 
(see ( [4.47D ) is obtained by the transformation ( [4.46| ). Using Af^g'^^d^: alone, there is 
no way to tell the difference between these two mass spectra. In contrast, our variable 
m!^2^'^'^ does not suffer from this ambiguity, and according to our results from Sec. iA 
the solution is always unique. 

The third advantage of the subsystem variable M^g'^'^^ is related to the expected preci- 



sion on the determination of the masses. As we pointed out in Sec. |4.4| and illustrated 
explicitly in Fig. |^, the kink A0(^'^'°) in the M^2^^lx(-^0' Pt ~ ^) function is much 
sharper than the corresponding kink AB^^'^'*^^ in the -^'r2'^da;(-^0)Pr = 0) function. 
This can also be seen explicitly from the two examples shown in Fig. |6|. As a result, we 
expect that the kink structure can be better identified in the case of M^2^^ax' which 
would lead to smaller errors on the mass determination. 
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('2 2 0) 

Of course, one could (and in fact should) use the experimental information from both Mj,^ max 
and m!^2^^1x^ available. Our main goal here is simply to point out the obvious advantages 
of the subsystem variable M!jp2^'^\ which so far has not been used in the literature. 

5.3 Hybrid method: Mt2 endpoints plus an invariant mass endpoint 



As discussed in Sec. 2A, any cascade with n > 2 will provide a certain number of measure- 
ments (2^) in addition to the M'T2 measurements discussed so far. In particular, for the 
n = 2 example considered here, there will be one measurement of the endpoint of the -Mr^^jj 
invariant mass distribution. The formula for the endpoint Mxix2,max in terms of the unknown 
physical masses Mq, Mi and M2 is in general given by 

1 



E^m = Mx,x,,max = Jf^^f {^1 " Mf) {Mf - M^) = {I - y)(l - z) . (5.16) 

In the case of tt events considered here, this is simply the endpoint of the invariant mass dis- 
tribution nifji of each lepton and its corresponding 6-jet. This distribution (unit-normalized) 
is shown in Fig. pi Unfortunately, here one is facing the same combinatorial problem as with 

(2 2 0) 

the Mrp2 ' variable - we cannot easily tell the charge of the 6-jet, therefore a priori it is not 
clear which 6-jet goes with which lepton. Fortunately, there are only two possibilities: the 
result from the correct (wrong) pairing is shown in Fig. |9| with the green (blue) dotted line. 
We see that the green histogram with the correct pairing has an endpoint at the expected 
location 




m 



w) {^w - O 



= \^ "^^^^^ ^ = 153.4 GeV, (5.17) 



m 



w 



with a relatively small tail due to the finite width effects. More importantly, the (blue) 



distribution from the wrong pairings is relatively smooth, and as a result the endpoint ( 5.17 ) 
is preserved in the experimentally observable (red) distribution, which includes all possible 
bi pairings. 



Now we can add the new measurement ( ^.16 ) to the previously discussed set of measure- 



ments (5^-5^). We obtain a total of five measurements for the three underlying parameters 



Mq, Ml and M2, therefore there exist two relations among the measurements. The first 



relation is already given by (5.12) and does not involve the invariant mass endpoint ( 5.16| ). 
The second relation is given by 

EL = • (5.18) 

We can now consider a hybrid method, which would make use of the invariant mass endpoint 
( 5.16 ), plus any two of the Mt2 measurements (^-5.7). In principle, one again needs to be 



careful and make sure that the three used measurements are independent. Fortunately, as 
seen from eqs. ( ^.12 , 5. IS), the invariant mass endpoint is independent from any pair of 



Mjn2 measurements. There are 6 possible pairs among the Mt2 measurements (^.4[-5.7), and 



in principle each one can be used in combination with the invariant mass endpoint ( 5.16| ). 



What is the best choice? We find that in all 6 of those cases one obtains a unique solution for 
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Figure 9: Unit-normalized mu invariant mass- squared distributions in dilepton tt events. The green 
(blue) dotted line corresponds to the correct (wrong) pairing of the leptons and the b-jets, while the 



red solid line is the average of those two distributions. The endpoint (5. It ) of the mu distribution is 
marked by the vertical red arrow. 



the masses Mq, Mi and M2. Therefore, the optimal choice is dictated by the experimental 
precision on each of the measurements (^.4-|5.7D. We expect that the measurement ([4.3|) of 



^T2^max ^^^^ ^^^^ precise due to the smaller cross-section for W~^W production. Similarly, 

max suffers from the combinatorial problem already mentioned earlier. Therefore for our 

illustration of the hybrid method we choose to use the m!^^''^^ endpoint the M^2'^'°^ 

endpoint ( jf)-7\ ), and the invariant mass endpoint (5.16). The solution for the masses in terms 
of those three measurements is given by 

i_ 

V2£^22i ( 2£'22i-E'iio + E22iEfm - Ef^y/ E'221 + 4-E'lio) 

Mo = ^ ^ - ^— , (5.19) 

^221 + "^Efm - E221 V E^2i + ^E^w 

\pl E22I Eim ( -E'221 \/ -E'221 + 4-^210 ~ -^221 ) 

Ml = ^ ^ ^ — ^— , (5.20) 

-^221 + 2-^2^ - -E'221 V EI21 + 4-E|io 

M2 = ^^^^^i%_^= . (5.21) 

-E'221 + 2-^^^ - -E22I V £^221 + 4-E'2io 

It is easy to check that substituting the measured values of the endpoints -E'221, -£-210 and 
Eim from (^), ( ^.11 ) and ( 5.17| ), into the three equations above yields the values for the 



neutrino, W and top quark mass, correspondingly. Notice that this determination was done 
in a completely model- independent fashion and without any prior assumptions, unlike the 
analysis of Ref. [40], which made use the known values of the W and/or neutrino masses. 
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Before concluding this section, we should emphasize again that the W+W and tt data 
required to test each of the three methods described in this section, already exists at the 
Tevatron. Therefore we use this opportunity to encourage the CDF and DO collaborations 
to perform a model-independent mass determination analysis along the lines presented here. 



6. Summary and conclusions 

In this section we shall briefly summarize our main results. 

• We compared the three main methods previously proposed for mass measurements in 
cascade decays with semi-invisibly decaying particles. We showed that the endpoint 



method (Sec. 2.1) and the polynomial method (Sec. p.2p are able to completely deter- 
mine the mass spectrum only if the decay chain is sufficiently long, namely n > 3. The 
same conclusion applies if we consider a hybrid combination of the endpoint and poly- 
nomial methods. As a corollary, when the decay chain happens to be relatively short, 
n < 2, these two methods are not sufficient, and one must somehow resort to the third. 



Mt2) method, in order to completely pin down the mass spectrum. Then in Sec. 2.3 
we argued that the Mt2 method hy itself is sufficient for a complete mass spectrum 
determination, even in the problematic cases of n = 1 or n = 2. In Sec. |5| we backed 
our claim with two explicit examples: W~^W~ pair production, which is an example of 
an n = 1 chain, and ti pair production, which is an example of an n = 2 chain. We 
showed that the Mt2 method in principle provides more than enough measurements 
for the unambiguous determination of the complete mass spectrum. 

When applying the Mt2 method, we generalized the concept of Mt2 by introducing var- 
ious subsystem (or subchain) m!^^^''^^ variables. The latter are defined similarly to the 
conventional Mt2 variable, but are labelled by three integers n, p, and c, whose meaning 
is as follows. The integer n labels the "grandparent" particle originally produced in the 
hard scattering and initiating the decay chain. We then apply the usual Mt2 concept 
to the subchain starting at the "parent" particle labelled by p and terminating at the 
"child" particle labelled by c. In general, the "child" particle does not have to be the 
very last (i.e. the missing) particle in the decay chain, just like the "parent" particle 
does not have to be the very first particle produced in the event. The introduction of 
the M^2^'^^ subchain variables greatly proliferates the number of available MT2-type 
measurements, and allows us to make full use of the power of the Mt2 concept. 

In Sec. ^ and Appendix |A| we provided analytical expressions for the endpoints of all 
n < 2 subsystem M^g'^' variables, as a function of the corresponding test mass Mc, 
and for arbitrary values of the ISR pT- Such results for general pT are being presented 
here for the first time. Only some special cases of our results have so far appeared in the 
literature, for example the functions M^2^''^^(Mo,Pt = 0) and m!^2'^'^\Mq,pt = 0) were 

{'2 2 1) 

already derived in [32] , from where one could also deduce the form of M^-g ' (^i > Pt = 



0). Our result for m!I?2^'^\Mq,pt = 0) discussed in Sec. 4.4 is also new. 
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In Sec. ^ we showed that the endpoint functions M^^''^\Mc, pt) niay exhibit up to 
three different types of kinks. All three kinks appear in the same location: Mc = Mc, 
at which point the value of the subchain endpoint Mj.2max coincides with the true 
parent mass: 

KTr^tiMc = M,,pt) = Mp . (6.1) 

(This equation generalizes eqs. (|4.7D , ( [4.38| ) and ( |4.53D .) However, the physical origin 
of each type of kink is different. 

1. The first type of kink was originally identified in [30,31] and arises solely as an ISR 
effect, so that in principle it should always be present at some level. The sharpness 
of this kink depends on the transverse momentum pj- of the ISR objects. This 
particular kink type was responsible for the kink feature AG*^^'^'*^) observed in our 



very first and simplest example in Sec. 4.1. There we generalized the existing 



analytical formulas for the endpoint M^2^^ix including the effects of the ISR. 
This in turn allowed us to analyze the amount of kink AG^^'^'") as a function of 
the ISR pt and the mass spectrum, see Fig. |5[ 

The second type of kink, encountered in Sec. |4.3| , was originally discovered in 
[29, 32] and is due to the variable mass of the composite system of visible (SM) 
particles used in the Mt2 calculation. To be more precise, this type of kink requires 
the following relation between the parent and child indices: 

p>c + l, (6.2) 

and does not require any ISR, i.e. it is present even when px = 0. Among the 
four M^g'^''^^ variables discussed in Sec. ^ only M^g'^'^^ satisfies this condition. 

(2 2 0) 

Not surprisingly, we encountered this kink during our discussion of the Mrp^^g^^ 



endpoint in Sec. 4.3. There we quantified the amount of kink AG^^'^'*^^ as a function 



of the mass spectrum, as shown in Fig. y(a). 



3. The third type of kink, which we encountered in Sec. 4.4 during our discussion 



of the M^2^'^^ variable, is new, and to the best of our knowledge has not been 
discussed in the existing literature^^ . This kink arises due to the decays of heavier 
particles above the parent level and exists even in the absence of any ISR: its 
presence simply requires the following relation between the grandparent and parent 
indices 

n > p , (6.3) 

while the child index c is arbitrary. Because of the upstream objects, this sit- 
uation does not correspond to the balanced momentum configuration discussed 



^■^Ref. [48] considered the decay chain q ~^ X2 ~* ~^ Xi ^-nd used the two leptons to form an Mt2 variable 
which in our notation would correspond to M^j'^''^' , i.e. n = 3, p = 2 and c = 0. Since this case satisfies 



both conditions (p^) and (6.3), the kink observed in [48] is a combination of the second and third kink types 



according to our classification. 
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in [28,32]. Therefore, the analytical expressions for A'It2 derived in those two pa- 
pers are not applicable and will need to be generalized [49]. In general, there are 
two sources of momentum imbalance in this case: upstream objects from grand- 



parent decays, as well as genuine ISR with px / 0. In Sec. 4.4 we concentrated 
on the former effect and provided analytical expressions for the endpoint func- 
tion M!j?2^^^^{MQ,pT = 0) and the associated kink A@^'^'^'^\ which allowed us 
to quantify the sharpness of the kink as a function of the mass spectrum, see 
Fig. |^(b). Comparing Fig. ^b) to Fig. |^(a), we saw that the new kink Ae(2'i.o) IS 
in general much more pronounced, and therefore offers better prospects for exper- 
imental detection. The corresponding formulas for the more general case, when 
both grandparent decay products and genuine ISR with arbitrary px are present, 
are listed in Appendix 

Of course, there are cases when two or even all three of these kinks will be simul- 
taneously present. For example, Mj,^ max^-^o^PT 0) will exhibit kinks 1 and 2, 

^r2,ml(^o,PT / 0) will exhibit kinks 1 and 3, while M^^a'.m"! (^o , ;5t / 0) will ex- 
hibit all three: 1, 2 and 3. 

Our Mt2 analysis in Sec. ^ revealed that in the case of an n < 2 cascade, there exist 8 
different measurements of subsystem Mt2 endpoints 

1. One measurement ( [4.5| ) of the endpoint M^p^^^^^ at zero test mass Mq. 



A zi ['2 2 1') 

2. One measurement ([4.25| ) of the endpoint Mj^^max 'z-ero test mass Mi. 

3. One measurement ( |4.43| ) of the endpoint M^g'^^i^ ^^'■^ mass Mq. 

^(2,2,0) 
'^T2,max 



4. Two measurements ( |4.44D and ( 4.45| ) of the endpoint Af^g'^^ix t'^o large values 



{Eb and E"^) of the test mass Mq. 

5. One measurement ([4.56| ) of the endpoint M^g'^^^^ at zero test mass Mq. 

6. Two measurements ( [4. 571) and ( 4.58| ) of the endpoint M^2^^lx large values 



{Ef, and E"^) of the test mass Mq. 



In addition, we also have one measurement of the endpoint ( 5.16| ) of the invariant mass 



distribution Mx-^x2i bringing the total number of available endpoint measurements to 
9. Given that an n = 2 chain contains only three unknown masses Mq, Mi and M2, 
it should be clear that the spectrum can be fully determined. In Sec. |^ we proposed 
three different methods for mass determinations, depending on what subsets of all 
those measurements one decides to use. For example, the pure Mt2 endpoint method 

r(n,P,c) 



of Sec. makes use of different Mj^^ maxi-^c = 0) endpoint measurements at zero test 



mass. The method of Sec. 5^ makes use of a single M^fi^^^axi-^^) function, measured 



at different values of Mc. Finally, the hybrid method of Sec. |5.3| combines some of the 

M^rdxiMc - 

endpoint M^j^^^^^ax- 



M^2'max(-^c = 0) cndpoint measurements at zero test mass with the invariant mass 
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In conclusion, our work shows that, at least as a matter of principle, the Mt2 concept, 
when properly generalized to include the subsystem variables M!jp^^''^\ can allow the mea- 
surement of the masses of all particles in SUSY-like events with arbitrary decay chains at 
hadron colliders. 
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A. Appendix: Analytical expressions for M^2^'^^^(Mc,Pt) 

The purpose of this Appendix is to collect in one place all relevant formulas for the various 
subsystem Mt2 endpoints M^P^^^lxi^c^Px) in the presence of initial state radiation (ISR) 

with arbitrary transverse momentum px- In all cases, we will find that M!j?2'^ax(^c, Pt) is 
given by two branches: 



K2ZLiMc,PT) = { \ , (A.l) 

Fjf'P'^^^Me.pr), if M,>M,. 

In what follows we shall list the analytic expressions for each branch pj^'^'^^ and F^'^''^\ for 
all possible {n,p,c) cases with n — c < 2. The grandparents Xn, the parents Xp and the 
children Xc are always assumed to be on-shell. However, any intermediate particles X^ with 
n>m>poTp>m>c may or may not be on-shell, and the two cases will have to be treated 
differently. Such an example is provided by the endpoint function m|^'^]J^^^(M„_2,pt) 



discussed below in Section [A.2| . For convenience, our results will be written in terms of the 
mass parameters ^(n,p,c) defined in ([4.1| ) 

Mn ( M2\ , 

These parameters represent certain combinations of the masses of the grandparents (M„), 
parents (Mp) and children (Mc), and do not contain any dependence on the ISR transverse 
momentum pT- As we discussed in Sees. ^ and ^ these are generally the quantities which 
are directly measured by experiment. Therefore, with the M'Y2 method, the goal of any 
experiment would be to perform a sufficient number of ^-parameter measurements and then 
from those to determine the particle masses themselves. 
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In some special cases, namely n = p, we shall also define pT-dependent /i parameters, 
where the pT dependence is explicitly shown as an argument: 



c){pt) 



1 + 



PT 



PT 



2MnJ 2Mn 



(A.3) 



When pt = 0, the py-dependent parameters (A.3) simply reduce to the p^-independent ones 
&■ 

IJ-{n,n,c){PT = 0) = M{n,n,c) • (A.4) 

We also remind the reader that test masses for the children are denoted with a tilde: Mc, 
while the true mass of any particle does not carry a tilde sign. 

A.l The subsystem variable M^!^'^"^ ^■*(M„_i,pt) 

The corresponding expressions were already given in eqs. ( 4.31| ) and ( ^.321 ) and we list them 
here for completeness: 



p{n,n,n—l) 



{Mn-UPT) 



F\ 



(n,n,7i— 1) 
R 



n,n,n-l) (Pt) + ^ (/^(n.n.n-l) (Pt) + + 



Pt 



(A.5) 



(M„_i,Pt) 



/^(n,n,n-l)(-Pr) + \j ( /^(n,n,n-l) ("Pt) - ^) + 

where the pr-dependent parameter /U(n,n,n-i)(PT) was already defined in (A.3): 

/^(n,n,n-l)(PT) = At(n,n,n-1) | ^ ~^ 



(A.6) 



PT 

2M„ 



PT 
2Mr, 



(A.7) 



As already mentioned in Sec. 4J, the left branch i^^"'"'" corresponds to the momen- 



tum configuration [p^'^ TT P^t^ TT Pt, while the right branch i?^"'"'" ^) corresponds to 
Puh^pS) Upt. 



A. 2 The subsystem variable M^g'^^"^ ^^(M„_2,pt) 

In this case there is an intermediate particle Xn-i between the parent X„ and the child Xn-2 
(see Figs. |l] and Our formulas below are written in such a way that they can be applied 
both in the case when the intermediate particle X^-i is on shell (M„ > Mn-i) and in the 
case when Xn-i is off-shell (M„_i > M„). 
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In both cases (off-shell or on-shell) we find that the left branch of M^g'^^j, ^^(Afn-2)PT) 
is given by 



n(n,n,?i-2) / X 
'l '{Mn-2,PT) 



PT 



A^(n,n,n— 2) (pt) + \ 



where the pr-dependent parameter ^J'(n,n,n-2){PT) was already defined in ( A. 3 ): 



fJ'{n,n,n-2){PT) = M(n,n,n-2) | 1 + (^^^ 



PT 
2Mn 



Pt 



(Ai 



(Ai 



The right branch F^'^'"' ig given by three different expressions, depending on the mass 
spectrum and the size of the ISR px' 



'{Mn-2,PT) 



' r^(n,n,n-2) , s -r ^ MI-MI_^ 

Fl ^(M„„2, -pt) , if > Af,._: ' 



(A.IO) 



{n,n,n—2) 



R,off 



(M„_2,Pt) , if < V and AMn,n-2{pT) < M^„_,^„,, 



l<r"'^(M„-2,PT), ifPT<^v?;;zi 



and AM„,„_2(pt) > M^„-ix„ 



Here AM„^„_2(pt) is a py-dependent mass parameter defined as 

2 



AM„,„_2(pt) 



M2 + ^ - M„_2 




which in the limit pT ^ reduces to 

AMn,n-2{PT = 0) = M„ - M„„2, 



(A.ll) 



(A.12) 



justifying its notation. Notice that AM„^„_2(pt) is always well-defined, since it is only used 
when the condition pT < (M^ — M^_2)/M„_2 is satisfied and the expression under the square 
root in ( [A.ll ) is nonnegative. The other mass parameter appearing in (|A.10|) , Mx„_-^xn,max, 
is the familiar endpoint of the invariant mass distribution of the {xn-i, Xn} SM particle pair: 



MtrJiM^ - M2_i)(M2_, - M2_2) , if M„„i < M„ , 

x„-ix„,max — ■\ [A.Lo) 
Mn - Mn-2 , if Mn-l > M„ . 



For example, in the special case of n = 2 and the intermediate particle Xi on-shell, eq. (|A.13| ) 

^{n,n,n-2) , j-,(n,n,n-2) 



reduces to eq. ( 5.16| ). The two expressions F^^'j and Fj^^^"' appearing in (|A.l[l|) are 
given by 



F, 



{n,n,n—2) 



R,off 



{Mn-2,PT) 



Pt 



Mn-2 + \ AM'{pT) + ^ 




(A.14) 
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{n,n,n—2) 



R,on 



{Mn-2,PT) 



^i_ix„,max +PlisiPT) + J + (pvisipr) - y) 



(A.15)' 

where AM„_„_2(pt) and Mri.^_-^xn,max were already defined in (|A.11| ) and ( A.13| ), correspond- 
ingly. The subscripts "off" and "on" in eqs. (|Al4| ) and (|A1^ ) can be understood as follows. 
When the intermediate particle Xn-i is off-shell and M„_i > M„, from ( [A. 11 ) and ( A. 15 ) 
we get 



Pt 



M, 



x„—iXn,max ■ 

(A.16) 



Now returning to the logic of eq. ( lA.lOl ), we see that in the off-shell case at low pT one 
would always use the expression Fji^^ff {Mn-2iPT) defined in eq. (lAll ), and never its 



alternative Pj^^" '^\Mn-2,PT) from eq. ( A.15| ). To put it another way, the expression 

^Ron^ '^\^n-2,PT) in eq. ( [A.15D is only relevant when the intermediate particle Xn-i is 
on-shell. 

Finally, the quantity PvisiPr) appearing in eq. ( A.15| ) is a shorthand notation for the 
total transverse momentum of the visible particles Xn and x„_i in each leg: 



In the case relevant for F: 



(n,n,n— 2) 



R,on 



Pvis — \PnT +P[n-1)T\ ' 

, the value of pms is given by 



Pt 

PvisiPr) = (/^(n,n,n-l) + l^{n,n~l,n~2)) l^JTr + l/^(n,n,n-l) " /^(n,n-l,n-2) I 1/1 + 



2M„ 



P\ 

AM?: 



(A.17) 



It is easy to check that in the limit of ^ our eqs. ( |A.8| ) and ( A.10| ) reduce to the known 
results for the case of no ISR (eqs. (70) and (74) in Ref. [32]). 



The left branch F-. 



{n,n,n—2) 



in ( |A.8| ) corresponds to the momentum configuration 



while the right branch in ( A.l[l| ) corresponds to 



In the latter case, F^^jr'^ "^^ is obtained when X„_2 is at rest: P^_2)t ~ ^' while F^'^^ 
corresponds to the case when P^_2)t ~ \pt ~ Pvisipr)- 

A. 3 The subsystem variable M!j?^'^^^"'~'^\Mn-2,PT) 



Here we generalize our pT = result (4.49) from Sec. [4.4| to the case of arbitrary ISR px- 



-,{n,n—l,n—2) 



(M„_2,Pt) 
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^(n— l,n— l.n— 2) {pt) + \ /^(n— l,n— l.n— 2) 



2 ^ ^ 

Pt 



(A.18) 



(n,n— l,n— 2) 



(M„_2,pt) 



^(„_l,„_l,„_2)(-j5Tj + ^/ ( ^(n-l,n~l,n-2)(-^r) " ^ ) + ^""2 



where we have introduced the shorthand notation 

PT=PT + 2/i(„,„,„_i)(pT) 



2 ^ 2 



(A.19) 



(A.20) 



Notice that the second term on the right-hand side contains the py-dependent /i parameter 
defined in 



The left branch F, 



(n,n— l,n— 2) 



in ( A.18| ) corresponds to the momentum configuration 



TT {piS TT pt) , 



while the right branch _p^"'" '^^ in ( A.19| ) corresponds to 



4l,)T n {pS TT pt) 



It is worth checking that our general pT-dependent results (|A.18 ) and ( A.19 ) reduce to 
our previous formulas ( [4.50 ) and ( 4.51 ) in the — > limit and in the special case of n = 2. 
First taking the limit px ^ from ( |A.20D and ( [A.7D we get 



lim = 2/i(„ „ „„i) , 

^li™p/^(n-l,n-l,n-2)(PT) = /^(n-l,n-l,n-2) (2Af(n,n.,n-l) ) = p(n,n,n-2) — /^(n,n.,n-l) 
^lim^ Ai(n-l,n-l,n-2)(— Pt) = ^(n-l,n-l,n-2) (" 2/X(„_„^„_i)) = H{n,n-l,n~2) ■ 



(A.21) 

(A.22) 
(A.23) 



Substituting ( [OIl - IOsD into ( [Al^ ) and ( [AloD , we get 



T(?i,n— l,n— 2) 



(M„_2,i3r = 0) 



/^(n,n,n-2) " ^(n,n,n-l) + y' ^^ln,n,n-2) + 



^^ln,n,n-l) 



(A.24) 



r-i(n,n— l,n-2) / „s. 



l^(n,n-l,n-2) + y (/^(n,n,n-l) " M(n,n-l,n-2)) + ^n- 



which are nothing but the generalizations of ( 4.5C| ) and ( 4.51| ) for arbitrary 



n. 
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